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Chapter 1 
Introduction 



Astrophysical observations show the existence of gravitational effect of an 
unknown type of matter that does not emitt electromagnetic radiation, thus 
being invisible to every direct relevation with the actual experimental instru- 
ments. For this reason this kind of matter has been called dark matter. 
There are many evidences of its existence: the rotational curves of the galax- 
ies differ from the theoretical predictions made taking only count of the visible 
matter, but they are explained if we add a dark matter contribution to the 
calculation. 

An estimate of the dark matter abundance in the Universe can be made 
studying the CMB (Cosmic Microwave Background). This is one of the 
goals of WMAP (Wilkinson Microwave Anisotropy Probe) experiment, that 
has given very important constraints on the Universe composition. Analyz- 
ing these data it has been calculated that the abundance of baryons and of 
matter in the actual universe are, respectively, Q^h'^ = 0.024 ± 0.001 and 
^Mh"^ = 0.14 ± 0.02, where QMh'^ is defined as Q.Mh'^ = The sum of 

Per it 

these contributions is about 20% of the Universe mass, so the existence of 
Dark Matter is necessary to justify these data. 

Many hypothesis have been made to try to explain the nature of dark mat- 
ter: primordial black holes, massive neutrinos, boson stars, brown dwarfs, 
non- luminous matter (black holes, non luminous gases, ...). However, one of 
the most natural explanation is given by the WIMPs and XWIMPs, particles 
appearing in some supersymmetric models. 

A WIMP (weakly interacting massive particle) is a massive particle that have 
only weak or gravitational interactions; an XWIMP (extra weak interacting 
massive particle) is a massive particle whose interactions are weaker than 
those of a WIMP, so are not of SM origin. 

Both types are not predicted by the SM. They appear if we extend the SM 
with supersymmetry. We will show that supersymmetry implies the exis- 
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tence of massive fermions called neutralinos, that are superpartner of W, Z, 
or Higgs bosons. We will show that a new symmetry, the R-parity, grants 
the stability of the lightest of these particles (the so-called LSP, lightest su- 
persymmetric particles) , that in this way becomes a candidate to explain the 
dark matter. 

We will start describing the Minimal Supersymmetric Standard Model (MSSM), 
focusing on certain aspects that will be useful for our purposes, as the neu- 
tralinos sector, that will be central in this thesis. Then we recap the results 
for the cross-sections of the annihilation XoXo ^ // and how these affect the 
LSP relic density calculation. 

Then we will examinate a different model, in which the MSSM is extended 
by an additional U{1) gauge group. We will describe the changes brought 
by this extension, especially in the neutralinos sector, then we will calculate 
the new interaction vertices and the related cross sections that affect the 
LSP relic density calculation. We will focus on the case in which the LSP 
is mainly a combination of particles introduced by the extension, while the 
LSP is mainly a MSSM neutralinos combination. Then we will modify the 
DarkSUSY package to perform calculations in this extended framework and 
we will discuss the results in the most general case. In the second part of 
the thesis we will study another phenomenological signature of our extended 
model, the asymmetry. This is one of the main observables that will be stud- 
ied at the LHC because can be used to impose constrains on the new physics 
beyond the SM and to distinguish among different models predictions. 
We will start describing the asymmetry at the LHC, where the simmetry 
of the initial state (pp) imposes cuts on the parameter space that lead to 
different asymmetry definitions. We will explain these definitions. 
We will calculate the cross section of the process pp — )■ e+e" and we will 
use these results to calculate the asymmetry at the LHC for all the defini- 
tions previously described. We will optimize this computation calculating 
the significance related to each definition and then choosing the cut on the 
parameter space in order to maximize its value. 

In this way there will remain only the dependence of the various asymmetries 
from the free charges of our extended model, so we will study its behaviour 
with respect to couple of them, showing bidimensional plot of the results. 
Finally we will study the general case, in which we will have dependence 
from three free charges. 
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Chapter 2 

MSSM dark matter 



2.1 MSSM 

2.1.1 Model setup 

In the SM the radiative corrections to the Higgs mass are proportional to 
the square of a renormahzation cutoff, that we will call A^y. This is an 
ultraviolet cutoff, needed to deal with the loop corrections of the theory. 
However, this implies that its value is large, much larger than the order of 
magnitude expected for the Higgs mass, that is around (150 GeV). This 
implies that the corrections are larger than the tree level value, making the 
Higgs mass larger and larger. Because the fermions and bosons masses in the 
SM are related to the Higgs mass, this problem will propagate also on their 
masses, affecting all the mass spectrum of the model. 

This problem can be solved extending the SM i.e. imaging that it is not valid 
at every energy scale but it is only the low energy approximation of a more 
general theory. One of the most promising extensions is supersymmetry. 
It is not a symmetry similar to those of the SM because it connects boson to 
fermion and viceversa. Thus the first request for a theory to be supersym- 
metric is that the number of fermions equals the number of bosons. 
The MSSM is the minimal possible supersymmetrization of the SM, i.e. us- 
ing the minimal number of particle to include the SM and to have a well 
defined theory. 

The first idea is to connect each boson to a fermion of the SM. Unfortunately 
the simple counting of the number of bosons and fermions in the SM shows 
that the first ones are 28: 
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while the fermions are 45: 
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Because supersymmetry associates a fermionic field to each bosonic field 
it is obvious that we cannot require this symmetry without changing the 
number of particles of the theory. Thus the MSSM is built imposing that 
each SM field is a component of a chiral superfields if it is a fermion or a 
bosonic superfield if it is a boson. In this way each particle has naturally its 
own superpartner. 

This method does not apply on the Higgs field: there is a general theorem 
that states that in a supersymmetric theory terms with ip"^ are forbidden in 
the superpotcntial. So it was introduced a second higgs with the quantum 
numbers of the complex conjugate of the first one. In this way we can combine 
the two higgses to obtain terms in the lagrangian that are invariant under 
the usual gauge symmetries. 

For completeness we give the higgs-fermion couplings after the introduction 
of the second Higgs field: 

W = XijU^HaUj + X',^U^H'eabDi + X'ljL'^H'eabEj + fiHaH'^ (2.1) 

where C/", U, D, L", E are the chiral superfields of quarks left, quarks up right, 
quarks down right, leptons left, electron (or mu or tau) right. Ha, are the 
two higgses. 

2.1.2 R-symmetry and R-parity 

The R-symmetry is a global C/(l) symmetry. Its action is defined in the 
sequent way: 

R(j){x, B) = e2"^"0(e-^"^, x) (2.2) 

in which 4>{x, 6) is a chiral supermultiplet, that we know contains a fermion, 
its bosonic superpartner and an auxiliary field. The action of the symmetry 
on the components is: 

A e^^'^A 

ijj e2^("-^)«-0 (2.3) 

P g2i(n-l)a p 
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where Ais the scalar component, ^jJ is the spinorial component and F is the 
auxihary field. 

The discrete version of this symmetry [1] is called R-parity and assigns a 
quantum number to each particle in the sequent way: 

= (_l)3(B-L)+2. (2.4) 

where B is the baryonic number of the particle, L is the lepton number, 
s is its spin. We can verify that for each SM particle and for the higgses 
Pr = 1, while for their superpartners Pr = — 1. If this quantum number 
is conserved, as happens in almost all supersymmetric models, there are 
interesting consequences in which we are interested and that will be one of 
the subjects of this thesis. The point is the sequent: the LSP, i.e. the lightest 
particle to have Pr = — 1, is absolutely stable. This is obvious: because Pr is 
conserved each particle with Pr = — 1 cannot decay in one or more particles 
with Pr = 1. So the lightest sparticle must be stable. 

Furthermore, if the LSP is electrically neutral, it interacts only by weak 
interactions with ordinary matter, becoming a perfect candidate to describe 
dark matter. 

2.1.3 Soft terms 

The soft terms are added by hands to the Lagrangian of the MSSM to break 
the supersymmetry without losing the propriety of absence of quadratic di- 
vergences and respecting gauge invariance. It was shown [2J that the soft 
terms permitted in the MSSM are: 

Csoft = -(^M,A"A'^ + ^a^^^Vi0,0fc + ^&''0.0,+tV. + /ix.) + (m2)!.0^>i (2.5) 
Z b Z 

where A" are the gauginos and ^^are the scalar fields. Note that the soft terms 
introduce masses for the gauginos, then there is a mass for the neutralinos, 
in which we are interested. 

2.2 Neutralinos 

Neutralinos are the superparters or the neutral bosons. Thus in the MSSM 
there are four neutralinos: two higgsinos and two gauginos. In the basis 
{B, W", Ha, H^) the neutralinos mass matrix is [T]: 
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(2.6) 
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gi e (72 are the coupling costants of the groups U (1) and SU (2) of the SM. Vu 
and Vd are the expectation values of the two higgses. Mi, M2 come from the 
gaugino mass soft terms in (12. 5p . fi comes from the higgsinos mass term in 
(12.11) . the terms proportional to gi, g2 comes from the Higgs-higgsino-gaugino 
couplings. 

The values of the masses and the related eigenstates are obtained diagonal- 
izing the mass matrix. So the LSP, xo, will be a linear combination of the 
four gauge eigenstates: 



The presence of these four components make possible many annihilation pos- 
sibilities for a couple of LSPs. Considering only the tree level there are the 
following channels: 



XoXo hh XoXo -> HH XoXo ^ hH 

XoXo AA XoXo hA xoXo ^ HA 

XoXo ^H+H~ xoXo^W^m XoXo ^ZH (2.8) 

XoXo Zh XoXo ZA xoXo ^ W^W~ 

XoXo ZZ XoXo ff 



h, H, A, H~^, H~ are the five higgs components that remain after giving mass 
to the gauge bosons (initially the higgses, each composed by two complex 
fields, consisted in eight degrees of freedom). H~^, H~ are the charged higgs, 
the others are neutral. A is parity odd, the others are parity even. 
In [3]we can find a complete calculation of the amplitudes for all the annihi- 
lation processes of the LSP at tree level. 



Xo = aB + bW^ + cHa + dH'' 



(2.7) 
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Chapter 3 



Minimal anomalous U(l) 
extension of MSSM 



String theory suggests [l]-[26] the existence of extra U(l) symmetries, so 
we want to extend the MSSM. We are going to introduce an extra U{1) 
symmetry. We will have an extra abelian vectorial multiplet and a Stiick- 
elberg multiplet, because in this model we use the Stiickelberg mechanism 
pS] to break the extra U{1) symmetry and because the Stiickelberg particle 
is needed to cancel the anomalies via the Generalized Chern-Simons mech- 
anism. So this extension is related to the fact that many string or GUT 
theories have extra symmetries with respect to the SM. If these symmetries 
appear in the high energy theory they can influence the low energy theory. 
Now we want to describe the features of the model that we will study in the 
rest of this thesis. A detailed treatment can be found in j27] . 

3.1 Model setup 

In this section, we discuss how to extend the Minimal Supersymmetric Stan- 
dard Model (MSSM) to accommodate an additional abelian vector multiplet 
V^(o). We assume that all the MSSM fields are charged under the additional 
vector multiplet V^'-^^ with charges that are given in Table ISTTI where Qi,Li 
are the left handed quarks and leptons respectively while t/f , D^, are the 
right handed up and down quarks and the electrically charged leptons. The 
superscript c stands for charge conjugation. The index i = 1, 2, 3 denotes the 
three different families. are the two Higgs scalars. 

Since our model is an extension of the MSSM, the gauge invariance of 
the superpotential, that contains the Yukawa couplings and a /x-term, put 
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Table 3.1: Charge assignment, 
constraints on the above charges 



Qe<^ 



-Qq - Qhu 
-Qq + Qhu 
-Ql + Qh^ 

-QHu 



(3.1) 



Thus, Qq, Ql and Qh^ are free parameters of the model. 
3.1.1 Anomalies 

As it is well known, the MSSM is anomaly free. All the anomalies that 
involve only the SU{3), SU{2) and U{1)y factors vanish identically. How- 
ever, triangle diagrams with t/(l)' current in the external legs in general are 
potentially anomalous. These anomalies ar^ 



f/(l)' - t/(l)' - f/(l)' 

f/(l)' - f/(l)y - f/(l)y 

U{1)' - SU{2) - SU{2) 
U{1)' - 5t/(3) - 5f/(3) 

u{iy-u{iy-u{i)y 



A^'^ = Y,Q 



f 

A^'^ = j:QfTr[TifT, 



(2)^(2)n 
k2 J 



(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 



"'^We are working in an effective field tlieory framework and we ignore troughout tlie 
paper all the gravitational effects. In particular, we do not consider the gravitational 
anomalies which, however, could be canceled by the Green-Schwarz mechanism. 
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where / runs over the fermions in Table IXTl Qf is the corresponding f/(l)' 



charge, Yf is the hypercharge and Tf 



(a) 

ka ' 



generators of the G^^) = SU{2) and G^^) 



2, 3; fca = 1, . . . , diniG^"'' are the 
: SU(?>) algebras respectively. In 



our notation TrfTj^^T^^"^] = \5jk- AH the remaining anomalies that involve 
f/(l)'s vanish identically due to group theoretical arguments (see Chapter 22 
of t29j). Using the charge constraints (13. ip we get 



3 [q^h^ + 2,Qh^QI + Ql- 3g|„ (Qi + 6Qq) } 



■ (3Qq + Ql) 



'QQh^ (SQq + Ql) 



(3.7) 
(3.8) 

(3.9) 

(3.10) 
(3.11) 



Notice that the mixed anomaly between the anomalous U{1) and the SU{3) 
nonabelian factors A^^'^ vanishes identically. 

Anomalous U(l)'s and the Stiickelberg mechanism 



We assume that the U{1)' is anomalous, i.e. f l3.7p -f l3TTT|) do not vanish. 
Consistency of the model is achieved by the contribution of a Stiickelberg 
field S and its appropriate couplings to the anomalous U {!)'. The Stiickelberg 
lagrangian reads 



2 



Y.bi^'^S Tr [W^-^W^'^^) + b^^^S W^'^ 

_a=0 

The Stiickelberg multiplet is chiral: 




S = s + iV2e^s + O^Fs - iea^'Od^s + -^O^Oa^d^^s - -O^e^Ds (3.13) 

We decouple the scalar component s in real and imaginary parts: s = a + i(f). 
We assume that the real part assume an expectation value from an high 
energy potential 

The S sector of the Lagrangian is: 
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(3.14) 



a=0 



, a=0 



1 



-0A(^V^a^A(°) + z(a)A(^V'^a^A(°) - ^F5A(i)aw 



sqrt2 



^PsX^D^''^ + ^y2A(^V'*a'^F(°V5 ) + (0 O 1) 



63 is a theory free parameter, 63"^ are parameters that have to be fixed in 
order to remove the anomaly. Note that 63 — because there isn't anomaly 

related to SU{3). 

Without details, we give the condition on the b'^^ to eliminate all the 
anomalies: 



1287r2 
b^% - ^^'^ 

3.1.2 Soft terms 

The soft terms of the model are: 



C 



'3847r2 
' 1287r2 



soft 



nMSSM 
*-soft 



I r<new 
"T '-soft 



The expression of the new introduced soft terms is: 



(3.15) 



(3.16) 



>CL7* = -^(MoA(°)A(°) + /i.e.) - ^(^V'5V'5 + h.c.) (3.17) 
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3.1.3 Kinetic mixing of the U(l)s 

In the equation 13.141 we can see the presence of mixing between the 2 super- 
fields of the U{l)s. We have: 



'A 4 2 ^ 



(3.18) 

with 6 = —Ab^2^ gQgi{a) . To diagonahze the kinetic terms we use the matrix: 

^ ^ f Cs \ f Vc 

where: 



)- \ -Ss 1 )\v, ' ^2-^9) 



Cs = (3.20) 

S'5 = 6Cs 



3.1.4 D-terms 



The D-terms of the model come from the kinetic terms of the chiral multiplets 
and from fl3.14p . Their Lagrangian is: 



a=0 a=0 



+2 



■ 2 ,(4) 
,a=0 



(3.21) 



a = 0, 1, 2, 3 runs over the gauge groups, Zi are the lower components of the 
z-th chiral multiplet (except that of the Stiickelberg field) and ka = 

1, . . . , dimG^"-* are the generators of the gauge group G^°'\ 
Using the equation of motion we obtain the expressions: 
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Dc 



/ / 

^2^f2hfgl {i,s [ClXc] + /i.e.] + 2^}!? g\ [^5 {Sl\c - Ss\b) + h.c\ 
^^hfg^gr \i)s {C,\b - ICsS.Xc) + h.c\ | (3.22) 

\WY. Yf\zf\' + 2V2b^^gl [i;s (Xb - SsXc) + h.c] + 



+V2hfgQgi [iJsCsXc + h.c] 
\ E {^.4(Tf );-^- + hfgl [y2^.Af ) + h.c\ } 

k 



(3.23) 
(3.24) 
(3.25) 



3.1.5 Higgs sector 

It is worth noting that in our model there is no stiickelber-higgs mixing. This 
is due to the fact that with our field content, = 1 supersymmetry and the 
gauge invariance with respect of the extra f/(l) there can not be a term that 
couple the stiickelberg field with the higgses (on the contrary to [31]). In 
fact the only permitted term that contains the Higgses and the Stiickelberg 
superfields is: 

W oc e~^^HuHd (3.26) 

Calculating the terms of the scalar potential generated by this superpotential, 
we can found that there appears only the real part of S, that in our model 
is fixed, while the imaginary part, 0, that is the propagating field, is not 
coupled with the higgses. After the electroweak symmetry breaking we have 
four gauge generators that are broken, so we have four longitudinal degrees 
of freedom. One of them is the axion, while the other three are the usual NG 
bosons coming from the Higgs sector. 

The Higgs scalar fields consist of two complex SU{2) ^-doublets, or eight 
real, scalar degrees of freedom. When the electroweak symmetry is broken, 
three of them are the would-be NG bosons G"^, G^. The remaining five Higgs 
scalar mass eigenstates consist of two CP-even neutral scalars h^ and one 
CP-odd neutral scalar and a charge -|-1 scalar as well as its charge 
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conjugate H with charge — 1|§ The gauge-eigenstate fields can be expressed 
in terms of the mass eigenstate fields as 



K 



1 

71 



K 



R 



'/3± 



Vu 
Vd 

G 
H 



(3.27) 
(3.28) 



where the orthogonal rotation matrices Ra, Rpo, R/s^ are the same as in [T] 
Acting with these matrices on the gauge eigenstate fields we obtain the di- 
agonal mass terms. Expanding around the minimum f l3.27p one finds that 
/3o = /3± = /3, and replacing the tilde parameters we obtain the masses 



m 



hO,HO 



Umlo+ ({goXsf + ]{gl + gl) ] v 




4 



{goXsr + \{gl + gl 



+4 { {goXsY + -{gl + gl: 



v^m% sin2(2/3) 



(3.29) 



(3.30) 
(3.31) 



and the mixing angles 
sin 2a 



sin 2/3 

tan2Q; 
tan 2/3 



m 



HO 



+ m 



2 2 



m 



AO 



+ [{g,Xsf + \{gl + gl)) 



{{goXsY + M + gl)) 



m 



AO 



(3.32) 



Notice that only the and masses get modified with respect to the 
MSSM, due to the additional anomalous U{1)' . 



3.1.6 Neutral Vectors 

There are two mass-sources for the gauge bosons: (i) the Stiickelberg mech- 
anism and (ii) the Higgs mechanism. In this extension of the MSSM, the 

^ We define = G+* and = Also, by convention, /i° is lighter than H° . 
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mass terms for the gauge fields are given by 



(2)m 



(3.33) 



C^, B^ are the lowest components of the vector multiplets V^, Vb- The gauge 
boson mass matrix is 



Ml 



i/l 4 



91-1 



(3.34) 



where = 16(7o&|C| + g^{%P')Xl and the lower dots denote the obvious 
terms under symmetrization. After diagonahzation, we obtain the eigenstates 



A.. 



Z' 



92B^^ + g^ 
V9l + 9l 



(3.35) 



9l 



+ 9oQhu 



^9l + 9l 



2„,2 



-C, + 0[glM-f,,] (3.36) 



a, 



(^g,B, - g^V^'^) + 0[gl M^^,] (3.37) 



and the corresponding masses 



Mt 



Ml 



Ml, 







(Q 



1 + 



7i + ^i) ^0^^ 



647r2M, 



V(o) 



(3.38) 
(3.39) 

v' + 0[glM-l,] (3.40) 



where My(o) = 463(70 is the mass parameter for the anomalous U{1) and 
it is assumed to be in the TeV range. Due to their complicated form, the 
eigenstates and eigenvalues of fl3.34p are expressed as power expansions 
in go and l/My(o) keeping only the leading terms. Higher terms are denoted 
hyO[glM-%]. 



The first eigenstate fl3.35p corresponds to the photon and it is exact to all 
orders. It slightly differs from the usual MSSM expression due to the kinetic 
mixing between V^'^^^ and V'^^\ 
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For the rest of the thesis, we neglect the kinetic mixing contribution since 
they are higher loop effects which go beyond our scope. Then the rotation 
matrix from the hypercharge to the photon basis, up to 0[gQ, Myf^-^] is 



Ofj. 



1 



V 
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92 

\f9\W2 



„ 9oQhu'" 
92 



\j9{+9i 
91 



\j9i+9i I 



(3.4i; 



where i, j = 0, 1, 2. 



3.1.7 Sfermions 

In general, the contributions to the sfermion masses are coming from (i) the 
D and F terms in the superpotential and (ii) the soft-terms. However, in our 
case, the contribution comes only from the Dq terms 

= { (csgoQH. + ISsg^^ (^^^) + ^^Csgohia) 1 1 {CsgoQf - Ssg^Yf) \yf 



where the y/ stand for all possible sfermions. 



3.2 Neutralinos sector 

We have already showed that, with respect to the MSSM, we have two new 
fields in the neutralinos sector: tps and X^^\ Subsequently we will have: 

^neutralino mass = -^(V'°)^M^V'° + h.c. (3.43) 
in which we have chosen the basis: 

i^) = {^s, Ac, Xb, A(^), hi hi) (3.44) 

Mj^f receives contributions from many terms of the Lagrangian: 

1) MSSM terms 

2) soft terms coming from C^ljf 



(3.42) 
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3) h — h — X' terms 

4) terms coming from Caxion 

5) D-terms 

So we have the symmetric mass matrix: 



/ Ms 

9 



M 



N 



tnsc 



\ 



^sB -^alb^? 











Ml 






M2 



-goVdXs goVuXs 



_9rvd 
2 

92Vd 
2 



2 

92"" u 







/ 

(3.45) 

The 4x4 matrix in the low-right corner is the mass matrix of the neutrahnos 
of the MSSM. The diagonal terms proportional to Ms and Mq come from the 
soft terms arisen with the anomalous extension. The terms with Cs and Ss 
comes from the U{l)^s kinetic mixing. The off-diagonal terms, except msc 
and msB come from the couplings h — h — X'. The Xs definition is: 



goXs — CsgoQn^ — -giSs 



(3.46) 



go is the coupling constant of the anomalous U{1), gi is the coupling constant 

of the standard U{1). 

The terms msc ^ ^sb receives many contributions from different terms of 
the Lagrangian. Their complete expressions are: 



msc = V2i 2 {Clglhf + Sjglh^^^ - CsSsg^ht^) {goXg Av^ + CgMy^o^a) 



+i [-2Ssgllf^'^ + Csgogihf) gi Av^ + ^My,, } (3.47) 

msB = ^ { {Csgogib^2^ - 2Ssgfb?^ {goXs Av^ + CsMy(o)a) + 6^'^^? A^;^} 

These expressions contain terms of order higher than 1 in the coupling con- 
stant. So they are negligible with respect to those in M^y that come from 
tree level. They can be neglected at the leading order. 
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Chapter 4 
LSP decay 



4.1 Pure stiickelino LSP annihilation to pho- 
tons 

We want to investigate the changes brought in the model by the anomalous 
extension. So we require that the LSP comes from the anomalous neutralino 
sector introduced adding the new U{1). We start requesting for simplicity 
that the 2x2 sector related to the extension is decoupled from the MSSM 
sector. 

4.1.1 Lagrangian and Feynman rules 

To get a pure stiickelino LSP we need some approximations. So we rewrite 
the neutralinos mass matrix at tree level in the anomalous U{1) extension: 



/ Ms. 

I 2 



N 



V 



m 

'5 



sc 



msB 











-miSs 





-goVdXs 


goVuXs 


Ml 





2 


givu 

2 




M2 


g2yd 

2 


2 









-fi 












(4.1) 



We want to find the limit in which the LSP is a pure stiickelino. We start 
searching the limit of decoupling between the anomalous and the MSSM 
sectors. Remembering formula fl3.20p the first observation that we made is 
that, at tree-level: 



5 ~ ^ ~ 1, ^5 ~ 0, ~ Qh^ 



(4.2) 
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Substituting in fl3.47p and eliminating higher order terms with respect to the 
tree-level, we obtain: 



msc = V2[2gob3] 
msB = 



(4.3) 
(4.4) 



To complete the decoupling remember fl4.2p . All the terms off the blocks 
2x2 and 4x4 different from zero are proportional to Qhu- 
So we impose Qhu = 0- Now the mass matrix has the form: 



M 



N 



/ Ms 
2 


2V2goh 














\ 




Mo 




















Ml 





2 


givu 

2 










M2 


92V d 
2 


g2Vu 

2 



















V •■■ 













/ 



(4.5) 



Furthermore we impose the LSP to be a pure stiickelino. Then Ms << 
Mq, that implies: M^sp ~ ^s- Now we are in the desired situation, so we 
can study the terms in the Lagrangian that contribute to the decay. The 
part of the Lagrangian that contains these terms is: 



C = -l4'^)[-!=tr(A('^)a^rF(»))^5 + h.c] 



(4.6) 



We want to express (14. 6 p in Dirac notation, so we introduce a convention for 
the gamma matrices: 



d'^ 




(4.7) 



Using this rule, calculating the traces and remembering the formula '^^'^^F^i, 
[l^ il^]dnBy {Bu is the gauge field related to F^^) we obtain: 



1^(2) 

^^?i2(A^)^75[7^7l(5;.S.)^5 + ^^72'(A')^[7^7l(9,W^J)*5 (4.8) 



V2' 



V2' 



The interaction vertex associated to this lagrangian is: 
C contains all the constants: 



(4.9) 
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C = ^gfJ^^wf a = 1,2 (4.10) 

Among the others, C contains 62°\ parameter coming from anomalous U{1) 
that imphes C << g(^a) [2Z]- The notation gi^a) and 62"'' means that the index 
can be both that of the hypercharge and of the group SU{2). w is a factor 
that, based on which group we are considering, will be: 

sinOw 

f is the gaugino coefficient in the mass eigenbasis. For a generical gaugino 
A(") it will be: 

A^'^) = c[^^Xi + c^\2 + c^\, + cf (4.12) 

Xi5X25X3!X4 are the MSSM eigenstates of neutralino mass matrix. / is the 
weight related to the gaugino that we are considering. 

The propagator, remembering that both stiickelino and gauginos are Majo- 
rana spinors, are: 

itT^P^ (4.13) 
Finally, the spin projectors related to Majorana spinors are: 

A± = (4.14) 
2m 
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4.1.2 Cross-section calculation 

The annihilation is represented in figure 14.11 The Stiickehno mass is Ms-, 





Figure 4.1: Annichilazione LSP-LSP 

while m is the mass of the generic exchanged gaugino. Applying Feynman 
rules to this diagram we have: 



u{p2)C'^^(y^''k2&ep{k2)- 



— m] 



(4.15) 



{Pi - kiY 

We define: = (pi — fci)^. Now, for generality, we calculate the product of 
two diagrams in which two different neutralinos are exchanged, obtaining: 



MM' 



tr\ 



24^2^/2 



(/f-m2)(/2-m'2) 



(4.16) 



Ms)l^a^''k2siUi - m)^^a^'ki^{^^ + Ms)k\a,,-i^iU i - m')klarsi^] 



Now we concentrate only on the trace evaluation. Using the identities: 



a'^ku = -{li,Y-k^^) 



-7m75 



(4.17) 
(4.18) 



in formula fl4.16p we obtain (we do not write the (1/2)^ coming from the first 
identity, assuming it semplifies the factor 2^ in ( I4.16P ): 



tr[{-i,, - Ms)iI^,Y - k'2)if, - m)iI^,Y - k\) 
X (-^, + M5)(^i7. - /^i.)(/i - ^')(^27p - k2p)] 



(4.19) 
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Before going on with this calculation we want to list the formulas that we 
are going to use: 



^^4^, = -24 (4.20) 

r4^7,^ = 4(^i?) (4.21) 

t'^^Wt. = -2{gi^4) (4.22) 

tr{4^) = 4:{AB) (4.23) 
tr{4^qLf) = A[{AB){CD) - {AC){BD) + {AD)BC)] (4.24) 

Using these relations we obtain: 

{^lY - kl){-i>i + MsM,!"" - h,) = -4^i(pi^i) (4.25) 
Substituting this result: 



We have used the relations: 

!i=i>i- ^2 (4.27) 

^2^2 = (4.28) 

^1^1 = (4.29) 

Let's evaluate the term with 8 gamma matrices: 

tr[{i>,li,Yl/>,Hil/>iH2lp] = (4-30) 

= 32{pih){p2k2){pik2) - 16{kik2){k2P2)M^s 
We now demonstrate that (p2^i) = (^1^2): 

Pi - ki = p2 - k2 ^ Pi + k2 = P2 + ki ^ (4.31) 
^pI + 2(pi/c2) + kl^pl + 2{p2ki) + kf^ (P1/C2) = {P2ki) (4.32) 

Calculating also the term with only 2 gamma matrices we obtain: 
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MM' = Q^C^C'\p,h) 

(M|-m2- 2(^2^:2)) (M|-m'2-2(p2)) ^ ^ 

[mm'(/c2fci)(p2A;2) + 2(pi/ci)(p2/c2)(pi/c2) - M|(p2fc2)(A;2fci)] 

To obtain the cross-section we have to multiply it with a suitable kinematical 
factor that relates AiM' to the cross section. According to ^32j the general 
relation between the amplitude and the cross-section is: 



da = ^-^dQ 

F = A{\^\EA + \pt\EB) (4.34) 



dQ={2Tlf5^^\pC+PD-PA-PB) 



{2TiY2Ec{2'Kf2ED 



Indices A,B refer to initial particles, indices C,D refer to to final particles. 
Now we concentrate on the center of mass frame. So we have these additional 
relations: 



Pa = -PB 

Pc = -Pd (4.35) 

s = Ec + Ed = Ea + Eb 

The last equation comes from the second and the fact that, in our case 
mc = rn£i. Using these relations we specify the f l4.34p to our case: 



F = 4\^\{Ea + Eb) = l^lv^ (4.36) 



Now, we define: 



W^V~s = Ec + Ej, = {ml + l^n^/^ + {ml + I^H^/^ (4.37) 

Using the f l4.35p we can calculate: 

dW 1 \ „_>, dW EcEd 

1^1 — + — ^ d\^\ = ^ , g (4.38) 



d\^\ '""'V^c Ed) mEc + Eo 

Now, before re- writing dQ, we remember the formula of spatial differential 
in polar coordinates: 
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d^Pc = \pcfd\p^\dn (4.39) 
Substituting these relations and using the spatial part of the we obtain: 



So the complete pref actor is: 

-^M- (4-41) 

If we call ui and UJ2 the energies of final photons, Ei and £^2 the energies 
of the initial stiickelinos, we have in the center of mass: 

Pi = -P2 (4.42) 

So the differential cross section is: 

S = , y MMj (4.43) 

2{hk2)^{uji + uj2f (4.44) 

(^2^:2) = E2UJ2 + (^M| - £;|)a;2 cos ^ (4.45) 



(p2A;i) = E2UJ1 - {^Ml - El)ui cose (4.46) 

The indices i,j vary from 1 to 4, because we have gauginos and each of 
them generate 2 vertices, with slight different constant factors. These are 
the expression in each case: 



Ci = V2^(i)4^^c^a^^^ mi = ^xi (4-47) 

C2 = V2^(^)4^^c^&^^^ ^2 = m^, (4.48) 

C3 = V^c/(^)6^^^s^a^^^ mg = (4.49) 

The parameter are the coefficient of the neutralinos ex- 
pressed with respect to the mass eigenstates. 
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4.2 Pure Stiickelino LSP and NLSP coanni- 
hilation 



As we can see from formula (14.341) . the cross section of the stiickehnos an- 
nihilation is proportional to C^C'^. These are the typical couplings of the 
extra U{1) that are at least one order smaller than the electro- weak couplings 
of the MSSM. So we can expect that the cross section will be 3-4 orders of 
magnitude smaller than the electro-weak counterparts. We will show that 
this is unacceptable if we want to obtain a relic density estimate in agreement 
with the experimental results. So we want to study the situation in which we 
have a NLSP coming from the MSSM with a mass comparable to that of the 
stiickelino. As we will show in a following chapter in which we will study the 
Boltzmann Equation, this situation can lead to a relic density that satisfies 
the WMAP constraints. 

4.2.1 Lagrangian and Feynman rules 

The possible coannihilating processes are: 



^5A ^ if°if°^5A ^ /J°/i°^5A ^ A°A°^5A ^ H+H- (4.54) 



/ and / can be any fermion anti-fermion pair of the SM, and W~ are 
the electro-weak charged boson, and are respectively the lighter and 
the heavier parity-even Higgs, and H~ are the charged Higgses and 
is the parity odd Higgs. In this section we want to obtain an estimate for 
the cross section, so we consider only the process \I'5A — )■ // because it is by 
far the main contribution to the total cross section. We will study the other 
processes when we will deal with the general case. 

The \1/A7 and \1/AZ° vertices are the same ones of the stiickelino annihilation, 
(14.91) . Because for this simplified calculation we consider the exchange of 7 
and we call Cz the costant related to the exchange and that of 
the photon exchange. From the SM we already know the vertices in which 
appear the pair //. The vertex ffj is: 



^5A ^ // 
^sA ^ W^W- 
*5A ^ h^h^ 



(4.51) 
(4.52) 
(4.53) 



(4.55) 
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In this expression e is the unit of electric charge, Q is the electric quantum 
number of the fcrmion anti-fermion pair considered. 
The vertex ffZo is: 

-Yigv-gAin) (4.56) 



2 sin 9w cos 9w 

e is, as in the previous formula, the unit of electric charge, 9w is the Weinberg 
angle, gv and qa are the vector and axial couplings, which depend on which 
fermion we are considering: 
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9v 






1/2 


1/2 


e, yU, T 


-l/2 + 2sm2^VK 


-1/2 


U, C, t 


l/2-4/3sm2^VK 


1/2 


(i, s, b 


-l/2 + 2/3sm2^vi^ 


-1/2 



Now we write the propagator. The photon one is 

k"^ is the momentum of the virtual photon. 
The Z° propagator instead, after fixing the gauge, is 

where fc^ represents the virtual momentum. 

4.2.2 Cross-section calculation 

he two possible diagrams that we consider in this estimation of the \E'5 — A 
cross-section of the process ^^^A — )■ // are: 




Mp2) fim) 

Figure 4.2: LSP-NLSP coannihilation 

In this section we perform the calculation of the cross- sect ion. 
The amplitude is: 

M = -ik^vi'-f^[-i^,'^u]u2[eqfC^—urfpVA (4.59) 

I 9Zq ^ I - / Zq Zn \ ^ 
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(4.57) 



(4.58) 



The square modulus is: 



\M\ 



T 

■J- n 



9z„ 



C 



^0 h.2 



(4.60) 



and are: 



t; = = s^^^r (4.61) 

where B^jo. is the result of the trace over the initial particles, and Ai^ 
are respectively the results of the vector and axial part of the trace over final 
particles. Now we are going to calculate separately the traces. We start from 

= tT[{-i^ - MsM-i. - + mA)(^7a - iM (4.62) 

to obtain this formula we have used the 75 property: 



757m = -7/.75 
(75)^ = 1 

Now we develop the product in the trace: 



(4.63) 
(4.64) 



The second term we want to calculate is vl^": 
Again we perform the product to obtain: 

The third term, A^ is: 

AT = -tr[{h + mf)Yl,{i>, - m/)757l 
After developing the products it becomes: 



(4.65) 



(4.66) 



(4.67) 



(4.68) 
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= trli^Y^l'"] + m)tr[Yl''] (4.69) 

Calculating the trace using the formulas fl4.23|) . fl4.24[) and performing the 
product, we obtain: 



+ 



-{PfPf) {PxiPf){PfPs) + {PxiPf-){Pf-Ps) 



+ m. 



{PxiPs){PfPf-) + 



-'^{P\iPf){PfPs) - iPXiPf)iPfPs) - {P\^Pf){PfPs) - 2{px,Pf)iPfPs) 



{P\iPf){PfPs) + {PxiPf){PfPs) m) - MiMs [mj - (p/P/)^] 

-{PfPf) {P\iPf){PfPs) + {PXiPf){PfPs) 



(4.70) 



Substituting these results in eq. fl4.60p . we have the complete expression 
of the amplitude. Now we multiply the amplitude for the prefactor (14.411) . 
In the center of mass we obtain: 



da 

In 



(4.71) 



The sum is over all SM fermions. is the color factor of each fermion and 
its value is 3 for the quarks, 1 for the leptons. 



4.3 Mixed LSP and NLSP cohannihilation 

In this section we want to study the general case, in which the LSP is not a 
pure stiickelino, but it is a mix of all neutralinos. We also want to deal with 
the mixing among neutral vectors and thus we will consider interactions 
that appear at the tree-level but we have not considered in the previous 
sections because they are numerically less important than those we already 
mentioned. 

4.3.1 Boson mixing and currents 

Inverting equations fl3.35p . fl3.36p and fl3.37p . we obtain: 

W'^" = M (4.72) 
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Defining an = Qq Qh, 



2v' 



we have: 



{ COs(9y 



M 



-9l 



l+an2(g2+gp 





l+an2 (32+5,2) 

cos{Ow) 
l+an2{g'f+g'p 

an {cos(9w)g2+sin(9w) gi) 
l+an2(g2+g2) 



l+an2(g2+g2) 



an cos{ew)y/gl+g. 
l+an2(g2+g2) 



l+an2(g2+g2) y 



(4.73) 

This leads to slight changes in the currents that we want to calculate. 
Now we use this mixing matrix to express the interaction lagrangian that 
involves bosons and their currents with respect to the mass eigenstates: 



-iguliMuA'' + M^,Z^ + Ph,Z'^) (4.74) 

-i52j2(M2iA^ + M22Z^ + M23^'^) 

-igojliM,,A^ + Ms2Z^ + M^.Z''^) 
Considering only the terms proportional to A'^ we have: 

-tig2M2^jl + g,Mujl)A'^ = 

= -le [^M,,3l + |Mn^ ) A^ (4.75) 

Experimental constraints on the neutral mixing imply that the photon cur- 
rent must be equal to that of the SM. This can be done easily, because the 
same constraints [33j permit us to not consider terms of powers greater than 
1 in an. Later we are going to show that this is all we need to keep the photon 
current equal to that of SM. For simplicity we perform the calculation using 
the electron, but it is identical for each SM fermion. Using ( I4.75p . we have: 

jT = 7^21 (-0 eil^^eL + ^Mn (-^6^7^61. - eR^.ei^ (4.76) 

where we have used the SM quantum numbers of cl and cr. Now we want 
to express left and right components with respect to the Dirac spinor, using 
the standard definitions: 

^L/B. = ^(1 T 75)^ ^ ^L/i?7M^L/iJ = *7/. (^^) * (4-77) 

Substituting these definitions in ( 14.76 p . we obtain: 

,|M„(-i.,Ji^)e-e-.„(i±-Vl (4.7S) 
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Now we can read the analogous of the vector and axial couplings of the SM 
and we can impose the desired equivalence: 

vec = - f--^M2i - -^Mn - -^Afn^ = -1 (4.79) 
2 V 2 e 2 e 2 e V ^ ' 

ax=- f--^M2i - -^Mn + -^Mn^ = (4.80) 
2 V 2 e 2 e 2 e V 

This impositions can be realized if and only if ^'2^21 = Qi^ii = c As 
promised, we now show that neglecting an powers higher than 1 is sufficient 
to achieve our goal: 

= g2M2i - giMn = 

I sin{9w) + an^ 9i\J gl + qD \ _ ( cos^Ow) + an^ ^2a/^i1^ + gl 
^\ l + an\gl + gl) ) \ l + an\gl + gl) 

^ g2sin{ew) - giCos{9w) g2sin{ew) = giCos{9w) (4.81) 
The latter is identical to the SM relation and also implies (substituting in 

In the anomalous extension the photon current is identical to that of the SM 
up to terms of second order in an, which are experimentally negligible. 
Thus the related Feynman rule will be the same of the SM: 




(4.83) 



Now we consider the terms proportional to Z^: 

-i{g2M22jl + giM^2jl + goM^2jt)Z^ = (4.84) 

7^ COs{ew)M22j^ + Mi2j^ + M32J 
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Again we try to reconstruct the SM structure of the current, so we have 

evidenced — Now we use the same calculations of the SM case to ex- 
cos (tiw) 

press the neutral current with respect to and j^™ (instead of jj), plus j° 
introduced by our extension. 

= cos{ew)M22jl + sin{9w)M,2jl + ^-^cos{ew)M,2f^ = 

92 

cos{ew)M223l + sm(evy)Mi2(j7 - jj) + -cos(^^)M32jJ = 

^ g2 

{cos{ew)M22 - stn{ew)Mu)jl + stniOw) Mi2jfr + -cos{ew)Mi05) 

92 

Now we want to eliminate the terms with powers higher than 2 in an. Using 
their definition in fl4.73p we note that there are simplifications in M12 and 
M22, while the M32 term remains unchanged: 

cos{9w)M22 - sin{9w)Mi2 = 1 (4.86) 
sin{ew)Mu = -sin'^iew) (4.87) 

Remembering that (j^'^)'^*^ = jj^ — sin'^{6y/)j'i^ this relations implies: 

jT = OD'"' + -cos{ew)M,23l (4.88) 

92 

where the latter term is the anomalous contribution. Now we have to cal- 
culate how this term changes the axial and vector couplings of SM. We only 
compute the variation, because obviously the analogous term in the SM will 
give the known contributions. Defining Ql, Qr as the anomalous U{1) quan- 
tum numbers of the left and right fermions respectively: 



—cos{ew)Mj,2 [QLeiliieL + QReRj^eR) 
92 



^cos{9w)M32 (^QLeJ^, (~Y^^ ^ ^ Qn^^l^i ^) " 

—cos{ew)M'i2 M 1 e7^e - I 1 67^756 1 (4.89) 



So we have: 



9v = + f^cos{ew)M,2Qv 

9A = 9f + f,cos{ew)M,2QA ^ ' 
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Where Qy^ Qa are definecO]: 

Qv = (4.91) 

Qa = (4.92) 
So the Feynman rule is the usuah 

7 




(4.93) 

It remains only to calculate the current of the Z'. The lagrangian term 

is: 

-t{g2M23jl + giM.^jl + goM,^f^)Z'^^ = (4.94) 

Again, neglecting the terms in fl4.73p of order higher than 2 in an, we obtain 
that Mi3 = an^ygf^^~g^sin{6w),M23 = —an^Jgl + glcosiOy/) and M33 = 1. 
As in the previous case we substitute with j^™ — jj^, obtaining: 

f,' = Jj + f ^M,3 - ^Mi3) jI + '-^M,,j- (4.95) 

V^o go J go 

For the quantum numbers this implies (omitting the calculation that is anal- 
ogous to the previous case): 

= Qv + T' (SM,3 - fM,s) + qeifM,, 

The related Feynman rule is: 

^Note that in p. II) are listed the right-handed charges of the anti-particles, so we have 
to change their signs to use them in these definitions 
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4.3.2 Neutralino-neutralino-boson anomalous vertices 

We have already studied the situation where the LSP is a pure stiickehno. 
However in the general case we can have mass eigenstates that are an unde- 
fined mix of gauge eigenstates. In this situation we have to account for this 
mixing in the interaction vertices. To achieve this goal we write the relevant 
terms of fl37H|) : 



1,(0) 

,{4) 

-^9li)^^'Yl^[r,mW^^s + ^9i9o[{X'^'Yl,b',r]d,C, +(4.98) 

(A(°))%[7^715/.5.]V'5 

This is the lagrangian in the gauge basis. Now we need the rotation matrices 
that maps bosons and neutralinos to the mass basis. For the bosons we have 
written the solution in the previous section, in formula (14. 73 p . For the neu- 
tralinos we cannot write an explicit matrix, because it is mo del- dependent. 
So we call the matrix and, ordering the basis as A^^-*, A*^^-*, A'^^, X'^^,-ips, 
we have: 

A('^ = 5^iVi>, = 5^iV*n, 

i i 

\^'^ = J2^>^ A(°) = J]iV>, (4.99) 

i i 

i's = '^N*^ni i = l,...,6 
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where rij are the mass eigenstates. Substituting this expression in (I4.99p we 
obtain: 

i j 
i j 

^^7o^?i 5^ rzf iV*,75[7^ 715^(^32^° + M,,Zl) N%n, + (4.100) 

1,(4) 

^^?o^7l ^^^6.75[7^ 715^(^11^- + ^12^° + M,,Zl) Y N*, 



This expression is very long, but if we collect the three different interactions 
it can be written in this more compact way: 



= 5ZC'A(^,J>f75[7^715M^.n, + (4.101) 

id 



where: 

CA{h3) = -^{bi'^ g!N*,N*,Mn + bf^ glN;,N;,M,, + 

-^gogiN;^N;,Mn + {i ^ j)) (4.102) 
Czit,j) = -^{bP glN;,N;,Ms2 + b^^^ glN*,N;,Mu + b^? glN^^N^.M,^ + 

1,(4) ,(4) 

-^g^g^Nl,N*,M^2 + -^gmN;,N*,M^2 + {i ^ j)) (4.103) 

^(0)^2Ar* Ar* ]\^_ , hW^'^AT* AT* A/T L ^(2)^2; 



^ib'^'g'oN;,N;,Mss + b\'> giNl^N*,M,, + b'^' giN*^N*,M,^ + 



.(4) ,(4) 
"2 „ „ AT* AT* n/r , ^^2 



-^g,g,N*,N;,M,s + -^gmN;,N*,M,, + (z o j)) (4.104) 

As we can see comparing fl4.102p with (14. 8 p we have the same structure of 
the lagrangian: the only changes are in the couplings just written. So the 
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vertices are all equal except that we have to use the general couplings. We'll 
have three possible interactions: 




(4.105) 

with respective Feynman rules: 

CAit,jh,[Y,YVK Cz{i.3)l,b^^l'']ik, Cp{^,3)^,Y1\Y]^K 

(4.106) 

in which /c^ is the momentum of the vector boson. 

4.3.3 Anomalous changes to higgs-higgs-boson and boson- 
boson-boson vertices 

We want to calculate the changes in the interactions of higgses and bosons 
caused by the anomalous extension. We start calculating the higgs-higgs- 
boson vertices. The lagrangian can be written in the same way of the fermion- 
fermion-boson one: 

(4.107) 

We have also an analogous term for Hd- Because the lagrangian is equal to 
that already studied, the currents will change in the known way when we use 
the mass eigenstates: 

jem _ -3 _|_ -y 

= 0-^)"'' + -cos{9^)M.. jI (4.108) 

92 

V^o 90 J ^ 90 
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In this expression we have the physical vector bosons, but the higgses are the 
gauge eigenstates. Using the formulas of [IJ, we obtain: 



HT 



cos(3 H+ 

-^{cosa + sina + -^cos(3 A° 

^ ^-sina h° + cosa + -^sin(3 
sin (3 



V2 



(4.109) 
(4.110) 



Remembering that (p d 



we know that if we have the same 



fields this derivative is 0. Using this fact we calculate: 

cos^[5 H~ d ^H^ + i cos a cos f3 d ^A^ + +i sen a cos 



H*, d = {H-y d ,H- + (if,")* d ,H', = 

^ y ^ y ^ y 

sin^P d ^H' - i sin a sin [3 d ^A^ + +i cos a sin (3 d ^^4° 

Using formulas fl4.109p and remembering that Qh^^ = 
action lagrangian becomes: 



hhv 



-I e A.H-'^nH' 



92 



cos {6 



w 



(4.111) 
(4.112) 

Qh^ the total inter- 



- - sm [iJw) 



^cos{ew)M32QH^ ] H-'d',H+ + 

92 



( - — —cos{9w)M32Qhu} {—sin a sinj3 — cos a cos f3)h^ d + 
V2 92 J 

cos{6w)M^2Qhu {—cos a sin[3 — sin a cos /3)ii° a (4.113) 

^ 92 J -I 



—i 


9oZ'^ 


+i 




+i 


(Qh 



1 /'5'2 



9i 



Qh„ + - -M23 - -M] 



2 \9o 



9o 



13 



9i 



^Mi3 ) H- d ,,H 
9o 



+i\QH - - ( — M23 - — Mi3 ) ) {cos a cos 13 + sin a sin (3)h^^ ^A^ + 
V 2 V^o 9o J J 



2 \9o 



^2^23 --M13 



9o 



[sin a cos (3 — cos a sin f3)H^ d ^A^ 



These are all the interaction terms between a vector boson and two higgses. 
The constant in each term is the constant of the respective vertex, whose 
Feynman diagram and rule are (from [5^): 

Now we are going to list all the couplings of these vertices. We call them 
F^'-^ , where A is the vector boson, BC is the couple of higgs entering the 
vertex: 
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k / 



Y 



(4.114) 



" cos{dw) 92 

V / V On On / On 



Ff^° = (4.115) 
^C^° = ^Hf^ - -M32QH^cos{ew)){-sm{a)stni(3) - cos(a)cos(/3)) 

Ff^" = 9o(Qh^ -l(-M23 - ^Mis))icos{a)cosm + szn{a)szn{P) 
V 2 V On On / / 



Ff^" = 



Fzo"^" = — rl-^fi - —M32QHuCos{9w)){cos{a)sin{f3) - sin{a)cos{/3)) 
" cos{6w) ^^92 ' 

F^"""" = 9o(Qh^ -l(-M23 - ^M,,))istnia)cosm - cos(a)szn(/3) 
V 2\go 9o ^ ^ 

Now we want to calculate the boson-boson-boson vertices. The lagrangian 
in the gauge eigenstates basis is the same of the SM: 



{^c^Vl - dpWDW^W^''] (4.116) 

As usual, to go to the mass eigenstates basis we have to use the formulas 
fl4.72p and (12. 6p . neglecting terms of order higher or equal to 2 in an. In this 
way the lagrangian becomes: 
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ie[{WlWp - WlW^)d''Af^ + {d^Wp - dpW^)W^^A'' - 
{d^Wl - dpWl)W^A''] + (4.117) 
ig2M2^[{WlWp - WlWa)d''Z'^ + {d^Wp - dpWo)W^^Z"' - 
{d^wl - dpWDW^Z'^] 

The first two terms give the vertices 7 — — W~ and Zq — — W~ and 
are exactly the same of the SM, so we obtain the well known Feynman rules: 




(4.118) 




(4.119) 



The third term, related to the Z', is due to the anomalous extension. 
However its structure is identical to that of the first two terms so its Feynman 
rule will be: 
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ig2M23[v"Hh - k2V + V^Hk2 - k^T + ^""(^3 - fci)^] 



(4.120) 

4.3.4 Calculation of coannihilations in fermions 

In this section we want to calculate all the extra U (l)-related contributions to 
the coannihilation of two neutralinos in a couple fermion-antifermion. The 
rules (l4.1U5p . f05D . fHl?^ and flOTD imply that there are three s-channel 
coannihilation processes: 




Summing this three contributions we have: 

M = -ik^vi^i^Yi ^,-iu\u2[eqfC A{id)-^u^1 pVi + (4.121) 
g2 rj^P _ rfP 

The tensorial structure is the same of fl4.60p so we have the same Ty and 
Ta- The difference with this latter expression is in the couplings and in the 
propagators, that gives us the amplitude square modulus: 




45 



Multiplying this result for the prefactor fl4.4ip we obtain the differential cross 
section in the center of mass frame: 



da 



1 


\pI\ 


647r2 


\Pi\s 



\M\^ (4.123) 



4.3.5 Calculation of coannihilation in vector bosons 

There are three diagrams for the cohannihilation in — W~ . They are 
analogous to that of the cohannihilation in / — /. The only difference is in 
the final states, which are a couple of vector instead of a couple of fermions: 




g2M23Cp{i,j)r] 
fc2 - M|, 



up 



-) [Vpaik - P3)[S + VafsiPs " Pi) p + VPpi.Pi " Pk)o\ 



a 13 



The square of this amplitude gives: 



lA/f|2_/iM"D r eCA(i,j) , g2Cos{9w)Cz{i,j) , g2M23Cp{i,j) ^2 . . 
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As for the fermionic final states, we have the product of two tensor terms. 
B^a belongs to the incoming neutralinos, and it is equal to fl4.62p . y4^° 
instead is one of the vectors in the final states. To calculate its expression 
we have to remember that, summing over the polarization p3], there are the 
identities: 

44 = Vap (4.126) 

Using this formula, and making the square modulus of the vector terms, we 
obtain: 



t^{k - psf + {k- psTiPs - p^r + (P4 - krik - p,r + (4.127) 

iP3 - PiTik - p;Y + 4(p3 - PATip-i - PaY + (P3 - PaT{Pa - kY + 
{k - PsYiPA - kY + (pa - kYiP3 - PaY + t'iVA - kf 

Computing the product A^^Baij, we obtain: 



|..|2 reC^ , g2Cos{9w)Cz g2M^23Cp^2 r , ^ , , 

-k^{2k^ {mxm^g + (P1P2)) + l{kp'i)mxm^g + 7{kpi)mxm^pg + 3{kp3){piP2) + 
3{kp4){piP2) + {kp2) {{piPs) + (piPa)) + {kpi) {8{kp2) + {P2P3) + {P2Pa)) + 
6(p3P4) {mxm^g + {P1P2)) - 8{piP2)m'l^ - I2m%,mxm^g + 4:{pip^){p2P^) - 
^{PiPa){P2Pz) - Q{P2P3){P2Pa) + '^{PiPa){P2Pa)) + {kpi){{p2P3) {k^ + ^kps) - 6{kp4)) + 
{P2Pa) {k^ - Qikp3) + 4(fcp4)) + 8{kp2) {{kps) + (kpi) - 2m^) + (4.128) 
12{kp2)ip3PA)) + k'^ikp3)mxm^^ + k'^{kpi)mxm^g + k'^ (kps) {piP2) + k'^{kp^)pi2 + 
k^{kp2){piP3) + {kp2){piPA) {k'^ - Gikps) + 4(/cp4)) - 2k^{pip2)m^ - 
2m^mxm^gk^ - 2{kp:ifmxm^g - 2{kpiYmxm^^ + Q{kp3){kp4)mxm^g - 
2{kp3Y{piP2) - 2{kp^fpi2 + Q{kp'i){kpi){piP2) + 4:{kp2) (kp^) (piPs) - 
Q{kp2){kp4){piP3))] 

As usual the differential cross-section in the CM frame is: 
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4.3.6 Calculation of coannihilation with final Higgs par- 
ticles 

Considering the coannihilation in a couple of Higgs fields there are different 
possibihties, with related diagrams: 




Where Hi and H2 can be the pair: H ,h A or H A . Using the 
Feynman rules fl4.106p and fl4.114p . we compute the amplitude: 



+ 



- Ml, 



+ 



(4.130) 



Now we want to calculate the tensor term, remembering the kinematic rela- 
tion = + p1: 



fc''t'i75[7^,7jM2?7''''(j93 +P4)p = 

Vi[Ji,li]ui = Q (4.131) 

So the cross-section is and there is not an higgs contribution to neutralino 
coannihilations. 
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Chapter 5 

Calculation of LSP relic density 



We want to give an estimate of the LSP relic density and then to compare 
it with the observed rehc density of Dark Matter. Now we want to describe 
the tools needed for the relic density calculations. We'll use a (+, — , — , — ) 
metric. 

5.1 Calculation method 

To calculate the relic density of a certain particle we have to consider the 
number of particles per unit of volume, n. We will also imagine to have coan- 
nihilation, that is the most general situation. This means that N particles 
are near in mass to the particle we are interested in. The evolution of the n 
in such a system is governed by the Boltzmann equation |35] : 

^ = -?>Hn - ^{(yijVij) {niTij - nf ) (5.1) 

aij are the annihilation cross section of the processes: 

XiX^ ^ // (5.2) 

the absolute values of relative velocity between i-th and j-th particle, 
defined by: Vij = \vi — Vj\ , rij is the number of particle per unit of volume 
of the i-th specie (index 1 refers to LSP), n^'^ are the number of particles per 
unit of volume of the i-th specie in thermal equilibrium, n = Ui, H is the 
Hubble constant and v is the relative velocity of the initial particles. 
The equation 15.11 is usually approximated considering rii/n = nl'^/Tf^. We 
also define the thermal average of effective cross section: 
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eg eg 



^effv) = Y,^a^,v,,)-^-^ (5.3) 



^eg ^eq 

Substituting we obtain: 

dn 



-3Hn-{a,ffv){n'-{n''^y) (5.4) 



{ceffVij) can be re-written in the sequent way [36] : 



i^effv) = . = (5.5) 



ntg Kg 



A is the annihilation rate per unit of volume at temperature T. We want 
to calculate it, using the classical gas approximation, that is the Maxwell- 
Boltzmann distribution: 



A = 5^(a.,.,>rnf = Y.W^J rfV.^^J>>-^'/^e-^^/^a,..,. (5.6) 

ij ij 

where Qi are the degrees of freedom of the i-th particle. Now we re-write A: 

^-l^J (2.)32E. (2.)32E, 

In this expression we have introduced: 

Wij = iEiEjaijVij (5.8) 

Now we define a function associated to the relative velocity, that is going to 
be useful in the calculation: 

V- E E- 

Pij = '^J- ^ Wij = ^Vij^soij (5.9) 

s is the Mandelstam variable that represents the energy in the CM frame, 
and it is defined as s = (pi +pj)^. Remembering that the relative velocity 
between two relativistic particles is: 



EiEj 

Substituting in (15. 9 p and re-writing with respect to s we obtain: 



= (5.10) 
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(5.11) 



2v^ 

From this relation and from (15.91) it is clear that Wij is a function of s only. 



Now we want to re-write the integral (15. 7p and change integration vari- 
ables to obtain a one-dimensional integral with respect to the variable s. We 
will use spherical variables for pi, modified spherical coordinates in which the 
z axis is generated hj pi, so the 9 angle is between pj and pi and the (pj angle 
is on the plane orthogonal to Pi. Re- writing the differential in this variables 
we obtain: 

d^Pid^Pj = \pi\'^d\pi\d cos 9id(l)i\pj\'^d\p'j\d(l)jd cos 9 (5-12) 

We can integrate all angular variables, except 9, because the integrand in the 
equation (15.71) is independent from them. So we have: 

d^Pid^Pj = 47r\pi\'^d\pi\4:7r\pj\'^d\pj\^dcos9 (5.13) 

Now we use the equation \pi\ = a/ Ef — mf and its analogous for \pj\ to 
change the variable in the differential from momentum to energy: 

d\pi\ = T^2EidEi (5.14) 

\Pi\ 

Using this result we obtain: 

d^Pid^Pj = 4:']T\pi\EidEi4:7r\pj\EjdEj^dcos9 (5.15) 

Now we want to change variables again. The target is to have a set of 
variables that includes s. So we define: 



E^ = Ei + Ej 

< E^ = Ei-Ej (5.16) 

s = mj + m'j + 2EiEj — 2\pi\\pj\cos9 

To express the differential in these new variables we have to calculate the 
Jacobian determinant: 
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dE, 


dE, 


dE+ 


dE^ 


ds 


dEj 


dEj 


dE, 


dE+ 


dE- 


ds 


d cos 9 


d cos d 


d cos 8 


dE+ 


dE. 


ds 



To calculate the determinant we have to invert the system fl5.16p . For the 
first two line the inversion is trivial: 



_ E^ (5.18) 



We clearly see that Ei ed Ej depend only on E+ ed E^. Given that \pi\ = 
Ef — mf, that \'Pj\ dependz only on ed if we invert the third 
equation of the system fl5.16p the only dependence from s is the explicit one: 

^s + mf + m] + l/2{El-E^_) 
cos 6' = r=rrr^, 5.19 

So we have = Now we can write down the known terms of the 

ds 2\pi\\pj\ 

Jacobian matrix: 

f 1 1 \ 
J=\ -2 -\ (5.20) 

\ d cos 9 d cos 9 —1 / 
\ dE+ dE. 2|p1||p5| / 

Using the third row we can calculate the determinant without knowing the 
remaining derivatives: 

detU) = —^^A -- + --) = (5.21) 

Using this result the equation fl5.15p reads: 

d^Pid^Pj = 271^ EiEjdE+dE_ds (5.22) 
The volume element is: 

d^Pi d^pj dE+dE_ds 



{2TiY2Ei {2TxY2Ej 8(27r)4 



(5.23) 



Now we have to find the integration limits. From the definition of spherical 
coordinates we know that —1 < cos 6' < 1. Now we want to find the condi- 
tions on Ei (the conditions on Ej will be analogous) starting from those on 

\Pi\- 
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\Pi\ >O^^Ef-mf>O^Ei>mi (5.24) 

Now, using Ei > rrii, Ej > rrij e \ cos9\ < 1 we have to find the conditions 

on E+, E- and s. 

The integration hmits become: 



s > {mi + mjY (5.25) 
> (5.26) 

E_ - E^^ ^1 < 2p,J (5.27) 

s \ s 

Now we can calculate A. The equation (15. 7p . in the new variables, reads: 

. f ..r e-^+/^dE+dE^ds 
A = J2 9,9^W.^ (5.28) 

ij 

We have already demonstrated that Wij is a function of s only, the exponen- 
tial is a function of E+ only, so we can easily integrate with respect to E^. 
It's sufficient using f l5.27p to derive: 



dE.=Apij\l^^^ (5.29) 



The integral becomes: 



. V- / ... El- se-^+/^dE+ds 
A = J2 9,9W.M^ V .4 (5.30) 



K2vr) 

These are the integral representation of the modified Bessel functions of the 
second type, that are crucial for the integration with respect to £'+: 



KJaz) = —— / dte-'H^-^^/W^ 



a? (5.31) 



We are interested in: 



dE^e-^-'^^El-s = T^sK, (5.32) 
£ dE,e-^^/^E,^E^-mj = (^) (5.33) 
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Using the first formula we perform the integration with respect to E^in 
(15.30 p . obtaining: 

^ = ^ E dsg,g,p,,W,,K, (^) (5.34) 



^3 

At this point we define Weff- 



9i9jPijWij = glPeffWeff (5.35) 

where: 



Peff = ys-4MI (5.36) 
Remember that Ms is the LSP mass. In other words, using (15. lip : 



^^^^ - ^ V s(s - AMI) ^ ' 

From equation (I5.36P follows ds = Speffdpeff, so: 

= / ^^^//^e//W^e//i^l (^) (5.38) 

Now we have to calculate the denominator of the equation (15.51) . using again 
the classical gas approximation and the relative Maxwell-Boltzmann statistic: 

-e. = E-r = E7|p/^'^>-"''' (5-39) 

i i 

Using the (I5.33P and (I5.14p . we go to spherical coordinates, obtaining: 

= E $ /" = 'i^K, {^) (5.40) 

Substituting (I5.40p and (I5.38P equations in (15. 5p . we obtain: 

r dpeffpLfWeffKi(^) 
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For example, if we have two particles coannihilating, the formula (15.371) reads: 




- AMI) ^ ^ ^ 

Indices 1 and 2 refer to LSP and NLSP. 

Now we have formulas for {a^ffv). We have to use this result to calculate 
the relic density. This is the normalization to 1 of the comoving number of 
particle per units of volume: 

Pcrit 

Pcrit is the density of comoving particles in universe. The expression for Vt^^ h? 
is calculated solving the equation (15.41) with appropriate approximations. 
The first step is to define Y = n/s, where s = S/B? is the entropy density. 
Assuming that the universe expansion is adiabatic, we have: sE? = cost, 
where R is the scale factor of the universe. We will also use the Hubble law: 
H = The dot means time derivative. 
We calculate: 

F = ^ - (5.44) 

s 

Furthermore: 



R = 3— (5.45) 



s J cost s 
Substituting this result in f l5.44p we obtain: 

Ys = h + 3Hn (5.46) 

So Ys is equal to the sum of the two terms in equation (15. 4p . that now, 
expressed in terms of Y, reads: 

Y = -s{a,fjv)iY'-Y,l) (5.47) 
We want to switch from the time derivative to the temperature one: 

dY _ dY dT _ dY dS dT _ dY rrn^dT 
'dt ~ 'dT~dt ~ drlldS ~ df^'^ ^dS 
In the last step we have used (I5.45p . The Boltzmann Equation now reads: 
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At this point we want to use adimensional variables, so we define: 

d^__M^ dY_ _dYM^ 

Substituting in equation f l5.49p we obtain: 
dY 1 m dS 



dx 3H x"^ dT 



a.ffv){Y'-Y') (5.51) 



We have to discuss the term -^j^- For this purpose we will use the Fried- 
mann equation for H and some parametrization for the functions S and p: 

H^=^-^ (5.52) 
3 

P = 9effiT)^T' (5.53) 

s = h,ff{T)^T' (5.54) 

geff and h^ff are effective degrees of freedon for the density and the entropy. 
They parametrize the deviation from the free gas behaviour multiplied by 
each effective degree of freedom. Using this relation we obtain: 



1 ds 1 I 3 30 I2n' fdKffiT). , 



3HdT 3 V SttG ^e//(T)7r2 45 V dT 



T-" + 3T'heff{T) 



^ heff{T) ^ T dheff{T) ^ 



AhG^^:^)'3Kjf{T) dT 
Now we define the effective degrees of freedom that we are going to use: 

1/2 _ heff{T) ( T dheffiT) 



^QeffiT) \3heff{T) dT J 

Experiments show that geff{T) and heff{T) are almost constant, so it is 

1 /2 1/2 

constant also g* , with a numerical value of gj ~ 9. However we will 
continue to consider those quantities as functions, dealing with the most 
general case. (We only omit to write the variables, so for example g^: = 
94T)). 

Using these results the Boltzmann equation (15.11) becomes: 
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Using this equation we can calculate the freeze-out point. It is defined as 
the point where A ~ SY^q, where 5 is a parameter of order 1 that is fitted 
to match experimental data. This fit gives the estimate 6 ~ 1.5. However 
we will use the symbolic expression for generality. The definition of A is: 
A = Y — Yeq. Now we rewrite the Boltzmann equation with respect to A. 
We have these identities: 



= ^ + ^ (5.58) 

dx dx dx 

Y^-Y^l = A{A + 2Yeq) (5.59) 



Substituting in fl5.57p we obtain: 



dA ^ dYeg / vrgf* 



Before the freeze-out Y ~ Y^g ^ A ~ 0, so the term with the derivative with 
respect to A can be omitted. Substituting the condition x/ ^ A ~ SYeq, we 
obtain: 



1 dYeg Ms Ing^ 



We start calculating the first member of the equation: for this purpose we 
write the explicit dependence from x of Y^g, using the equations (I5.40p and 
(1531: 



1 dY,g ^ _l__d_ /neg\ ^J_^45g^rmi_y_d_ I ^'-^2 (f^x) 

Y,g dx Yegdx \ S ) Yeg^ 47r4 \Ms) dx\ heff (^) 

(5.62) 

Now we calculate the derivative, using a property of the modified Bessel 
functions of the second kind: 

^^ = -i^._i(;.)-^ir.(^) (5.63) 
Using this equation we obtain: 
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Yeq dx Y,,^4n^Kjf[Ms) Kff {^) 
Substituting this result in equation f l5.6ip we have: 



Ms^l^Weffv)6{6 + 2)Y,, (5.65) 

In the non-relativistic case, Xf » 3 we can use an asymptotical expansion 
for the modified Bessel functions of the second kind: 

^.Woc + (5.66) 



In this limit this equation implies that Ki and K2 (^jj^^ are equal, 

SO; — ) nf \ — 1- If we assume that comoving entropy density is conserved 
we also have -^log {hefj {^)) = 0. The equation now reads: 



Ms^^{a,ffv)6{6 + 2)Y,, (5.67) 

X/ is the solution of this equation, that is numerically determined. For weak 
interactions the result is almost independent from effective cross section, and 
is Xf ~ 25. 

Now we can calculate the relic density of LSP. We are interested in the 
situation prior the freezing-out, when we are not near the equilibrium: this 
means Y » Y^q. Neglecting Yeq the Boltzmann equation becomes: 
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dY Ms I Tig* I /p-fioN 

'(TeffV)[Y ) (5.68) 



dx x2 V 45(7 

We can find the solution integrating from Tf to tfie actual temperature Tq 



Tliis is tlie complete approximate solution of the Boltzmann equation. 

To have an estimate we can assume that g* is constant and we can also 

neglect F/, obtaining: 

Substituting this result in f l5.43p we can calculate the relic density, remember- 
ing that (for our LSP) the relic density is given by p^g = MsSqYq, where sq 
represents the entropy density at the actual time. So we have the expression: 

n^,h' =P^ = (5.71) 

Peril Pcrit 

Using now pcrit = |^ and Mpi = together with (15.701) . we can estimate 
the relic density. There is also a very naive rule to estimate the relic density. 
It consists in the approximation of (cTe//f) as a constant; in this way the 
integration fl5.70p is trivial. The resulting approximate expression for the 
relic density is: 

^ , 2 3 X 10-27cm3s-i 

^^^^ ^ — ^r^A — ^^-^2) 



5.2 Relic density 

Our starting points are the annihilation and coannihilations cross- section al- 
ready calculated. If we consider only the annihilation, from equation fl4.34p 
we infer that the cross-section is proportional to C^C'^. Naive estimates 
from (I4.10p lead to a C and C of order of magnitude smaller than the typ- 
ical electro-weak couplings. This implies that the annihilation cross-section 
is roughly of the order of 10~^^cm^s~^. 

Substituting this value in (I5.72p we obtain VL^^h? ^ 10^ which exceeds the 
WMAP [27j constraints. We can exclude that the annihilation of an LSP 
coming from our extension could match the experimental data. 
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The case in which there are coannihilations is more interesting: the equation 
fl4.60p imphes that this cross section is proportional to C^e^. Thus, remem- 
bering that the masses of the LSP and NLSP are free parameters, we want 
to study if the LSP rehc density can be described by this process. 

5.2.1 Cross-section and Mandelstam variables 

The formula fl5.4ip implies that the calculation of the thermal cross-section 
requires that the total cross-section was expressed in Mandelstam variables. 
Their definition, with reference to the figure 14.21 is: 

s = {pi+ P2Y = {P3 + Pif = k'^ 

t = {p-i-pif = {pA-P2f (5.73) 

U = {Pi-Plf = {Pz-P2f 

Now we give some important properties of these quantities in the CM frame, 
defined by: 

Pi = -P2 (5.74) 
The momentum conservation implies: 

P3 = -P4 (5.75) 

So we have: 

{Pi + P2) = {El + ^2, 0) ^ s = (^1 + E^f = EIm = (Es + E^f (5.76) 

The explicit writing of the variables permits to calculate many other useful 
relations: 

s=pI+pI + 2(piP2) = Ml + ml + 2{pip2) = 2m) + 2(^3^4) 
t=pi+pl- 2(piP3) = M| + mj - 2{pip^) =ml + m)- 2(^2^4) (5.77) 
u=p\+p\- 2{pipi) = Ms + mj - 2{pip4) =ml + mj - 2(^2^3) 

These relations express the Mandelstam variables as functions of the masses 
and scalar products of the momentum of the external particles. In the scalar 
products there are the dependencies from the angular variables, that have to 
be made explicit in order to write the total cross-section from the differential 
cross-section: 
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= EiE2 + \pi 


2 




{pm) 


= EsEi + \p3 


2 




(pm) 


= E^Es - \pi 


\P3 


COS 6 


(pm) 


= EiE4+ \pi 


bs 


cos 9 


{P2P3) 


= E2E3 + \pi 


\P3 


cos 6 


{P2P4) 


= - \pi 


\P3 


cos 6 



In these formulas 6 is defined as the angle between pi and p^. At this point 
we also want relations that express all non-angular variables as functions of 
s. We want to solve the general case, in which the two initial particles (or 
analogously the two final particles) have different masses. We have to solve 
the system: 



ml + ml + 2E1E2 + 2\pi\ 



ml = Ef — 

"^2 = ^1 - IpiI^ 



The second and third equations give the solutions for Ei and £'2 with respect 
to that remains the only real unknown quantity. Solving the system 

we obtain: 

,r, s ml + m^ [ml — miY 
Now we want to calculate how this expression reduces if mi = m2: 



|pi|2 = |-m^ (5.80) 
^2 = El = s/A (5.81) 

Using these relations we can calculate the total cross-section and express it 
function of s. 

5.2.2 Relic density examples 

Now we want to study the relic density in simple cases of coannihilation. The 
simplest possibility is that in which the LSP is a pure stiickelino, coannihi- 
lating with a MSSM bino or a MSSM wino in a couple / — /. In this case 

is given by fl4.10p . We write 
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down the related total cross-section, using the formulas (14 . 71 p . f l3T^ . f l5 . 7lJ|) 
and integrating over the solid angle: 



- 2Mf + (4M| + s) Ml - QMssMi - 2M| + + M|s 



X 



127r (M|^ - s) s^/^^Mf - 2 (M| + s) + (M| - s)' 

(2mj + s) ( 2 cos 9 weqf{Ml^ - s) + sinOwgzoVfs) + {sin OwgzoCi'f)'^ {s 



This result can be used for a numerical calculation of the relic density. From 
the formula (15.51) . written in the case Ncoan = 2, we obtain: 

eq eq eq eq 

^,,^,^,,^ (''nv)/W2.v) + 2{.n^^^^ (5.83) 

where Q = 71^2 ■ The first term in the numerator can be neglected because 
the Stiickelino annihilation cross section is suppressed by a factor {C^°^ Y with 
respect to the MSSM neutralino annihilations (as already said) and thus 
(diif) ^ (o"22'i^)- The second term involves the coannihilation cross section 
and it is the termal average of (I5.83p . As we have seen in the previous section 
the anomalous couplings C.^ and Czq are both proportional to C*^^^ ~ b^2^9i, 
with 62^'* given by (I3.15P : 

where 63 = Mz'/Ago. With the assumption Mz' = 1 TeV as in [27j we can 
give the estimate: 

-q- ~ 5.76 X 10-'\3goQQ + goQif GeV-^ (5.85) 
Substituting this result in (I5.83p . given that the perturbative requirement is: 

gl ■ {^Qq + QlT < 16 (5.86) 

and given also that the typical weak cross section is {o'22v) — 10^^ GeV~^ 
we obtain: 
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in the case of a pure bino, while 



(q"i2f) 

{cr22v) 



< 10- 



5,10- 



(5.87) 



(5.88) 



in the case of a pure wino. So we can use the formulas fl5.83p and (15.721) to 
write: 



(2) 



l + Q 
Q 



n 2 



(1) 



(5.89) 



in which the relic density in the presence of coannihilations is related to that 
of the MSSM system without the stiickelino. We have calculated the MSSM 
relic density with the DarkSUSY package [SB] and then we have used it to 
calculate the coannihilating relic density for different mass gaps between LSP 
and NLSP. In figure (15. ip we have plotted the results. The red points satisfy 
the WMAP [37] constraints: 

0.0913 < VLh'^ < 0.1285 (5.90) 



Now we want to study the case = 3. This means that there are two 
particles with masses comparable to the LSP mass. This can be achieved 
very simply if the NLSP is a nearly pure wino: in fact the wino is almost 
degenerate in mass with the lightest chargino. Expanding all the terms in 
the sum (15. 5 p we get: 



cr. 



v) 



{cruv)jl + (cri2t;)7i72 + ((Ti3t;)7i73 + 
(0-21^^)7271 + (cr22f )72 + (0-23^)7273 + 
(0-31^^)7371 + (0^32^^)7372 + (0-33^)73 
(o-iit;)(nf )^ + (cri2t;)nf ^2'' + (o-i3t;)nf + 

{(y2iv)n^'rf^ + {<y22v){ri'^Y + {(^2i,v)n^n^ + 
{(y-iiv)nfnf + {(y^2v)nfnl^ + (0-331;) (ng''^^ 



(5.91) 



where in the last line we have neglected the terms (anf), (cri2f) and (cri3f) 
since these are the thermal averaged cross sections which involve the Stiick- 
elino. In the equilibrium approximation we can write: 



n. 



eq 



9i 



;i + A,)3/2e--/^» 



for i = 2,3 



(5.92) 
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Figure 5.1: Stiickelino relic density in the case in which the NLSP is a hnear 
combination bino-higgsino. Red (darker) points denote models which satisfy 
WMAP data. Left panel: A2 = 1%. Right panel: A2 = 5%. 



where gi are the internal degrees of freedom of the particle species and Aj = 
{rrii — mi) /mi. That permits to define and approximate the variables: 

eg 

Q^ = %i = -i^ + A,)3/2e--/^» for ^ = 2, 3 (5.93) 
^1 9i 

As in the N = 2 case, we use this variables to rewrite the termal cross-section: 



[1 + Q2 + Q 



3 J 
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This expression can be further simphfied assuming (7713 — m-zj/mi <^ ^/^fi 
QzlQ2 ^ ^3/52: 



1 + 



— 2('^MSSM^') (5.95) 



i + (i + S)g 



where ^ 

{(JUSSMv) = (0-22^') + 2 — ((723^) + ( — ) (cr33t') (5.96) 

92 \92/ 

Now we have to express this quantity with respect to a 2 particles thermal 
cross-section, in order to compute a rescaling factor as in the N — 2 case. 
The 2 particles cross-section in our case is given by: 

M = 

Mjnir + 2(a23^)n^V + {assv){niy 

Mjnl'f + 2(a23^)n>^'' + Mjnff 

_ {cr22v) + 2((723^)<?23 + (o"33^)(?i3 
(1 + ^23)^ 

3/2 



(5.97) 



where 



<523 = ^3/^2 = — U + 

92 \ m2 J 



j.^ "'3-"'2 
J mo 



(5.98) 



93 



~ 92 

since {m^ — m2)/mi ^ 1/^/ ^"^^ "^2 > ""^i then (ma — m2)/m2 ^ 1/^/- 
Note that the values of rf2-i ^'3 ^'^^ ^"^^ ^'^^ different with respect to those 
in the former case since now there are only two species in the thermal bath. 
We then find 



(<722^)+2|(a23t;)+ g) (a33t;) 
W^) - 7 72^ 

(cMSSMf) 



(i + l)' 



(5.99) 
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and inserting back this relation into fl5.95p we obtain 

2 



(3) \ 



1 + 23 
92 



Q2 



(5.100) 



The rescahng factor between the three and two particle species relic density 
is given by the following relation: 



2\(3) 



i + (i + S 



Q2 



2n(2) 



(5.101) 



In the same way of the N = 2 case we numerically calculate (i7/i^)*^^) with the 
DarkSUSY package, rescale the result according to (15.1011) and compare it 
with the WMAP [37] constraints for different mass gaps. We give the results 
in figures (15. 2p . 
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Figure 5.2: Stiickelino relic density in the case in which the NLSP is a wino 
while the NNLSP is the lightest chargino. Red (darker) points denote models 
which satisfy WMAP data. Upper left panel: A2 = 20%. Upper right panel: 
A2 = 10%. Lower left panel: A2 = 5%. Lower right panel: A2 = 1%. 
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Chapter 6 
DarkSUSY 



In this chapter we want to describe how the routines of the DarkSUSY pack- 
age permit to numerically calculate the relic density of the LSP and thus of 
the supersymmetric candidate to explain the dark matter abundance. Our 
aim is to use this package to calculate the relic density in the anomalous 
model. 

First of all we will describe how the routines of the standard package work, 
then we will describe how this routines can be modified to perform not only 
calculations in the MSSM but in our extended model. Finally, we will hst 
our results obtained with the extended package. 

6.1 Relic density with DarkSUSY package 

DarkSUSYis a program package for supersymmetric dark matter calcula- 
tions. It contains many set of routines, each with its purpose; the idea of the 
package is that the user can choose what set of routines is more useful for his 
calculation and write accordingly his own main program. We are interested 
in dark matter relic density and now we want to explain how the related 
routines work. 

6.1.1 Model defining routines 

• Common blocks 
In DarkSUSY the common blocks, which contain the variables that are 
used frequently in the routines, are stored in the folder " /include", 
that is in the DarkSUSY root directory. There are many files, all with 
the appropriated extension .h, that contain the common variables re- 
lated to a certain use. 
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The variables that are used to define the supersymmetric model are in 
the file "dssusy.h" . This file does not contain directly the variables, but 
calls four other files, "dsge.h","dsio.h", "dsdirver.h"and "dsmssm.h". 
"dsge.h" contains general numerical values, as tt, "dsdirver.h" contains 
values used by the program during the installation and are not used to 
perform calculations. Even "dsio.h"is not related directly to calcula- 
tions but sets the parameters that fix the unit in which certain routines 
will print their results, "dsmssm.h" is the most important file from the 
point of view of the calculations, because it contains the declarations 
and definitions of the numerical parameters that identify the particles 
in the routines, the definition of all the physical constants, such as cou- 
plings, Weinberg angles, mixing matrices, masses and so on. 
The variables that are used in the calculations of the cross sections are 
contained in "dsandwcom.h". The variables that set the parameters 
that will be used in the numerical relic density calculations are defined 
in "dsrdcom.h"and "dsrncom.h". 

• Initializing the program 

The program is initialized by the routine dsinit, that is located in the 
folder " /src/ini". This should be the first routine called in any main 
program, because it gives values to all the physical parameters such 
as particles degrees of freedom, quantum numbers, SM known masses, 
constants. Furthermore, in this routine the program sets the constants 
that decide how it will perform certain calculation. To explain this 
point, it is important to known that DarkSUSY contains many possi- 
bility to perform the same calculation, some quicker, others with higher 
precision. The user can choose the appropriate method for his needs 
simply changing the related value in dsinit. Furthermore the user can 
choose if the program has to consider loop effect and other details. 
Last, this routine call other subroutines that set constants used for 
calculations different from the relic density. We are not interested in 
those calculations so we do not deal with these subroutines. 

• Setting model-dependent parameters 

There are different ways to define a certain model in DarkSUSY. 
The first one is to introduce each single parameter of the new model. 
In this way the user writes directly the value desired for each parameter. 
This input is read by DarkSUSY in the " /src/su/dsgive_model.f' routine 
if we want an effective model defined at low-energy scale, or by " /src/rge/ 
dsgive_modeLisasugra.f'if we want an mSUGRA model whose low- 
energy parameters are calculated by RGE evolution. In the following 
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we will use the first option since we are not interested in the evolution 
from a high energy scale. 

The second way to define a model in DarkSUSY is to generate it ran- 
domly. The related routine is "random_model()"and can be found in 
the file " /test/dstest.f. This routine simply generates random values 
for all the model parameters. The range in which these parameters are 
generated can be manually changed by the user opening the routine 
and changing the desired minimum and maximum values. 
The last way to define a model is via an SLHA file. In this way the 
program reads the file using the routine " /src/slha/ dsSLHAread.f 'that 
sets the parameters according to the given file. 

• Defining all model-dependent quantities from the inputs 

There are three possibilities: if the inputs were given (manually or ran- 
domly) at high energy and calculated at low energy by RGE you have 
to use " /su/rge/dssusy_isasugra.f' ; if the inputs were given at low en- 
ergy (manually or randomly) you have to use " /src/su/dssusy.f'and if 
the inputs were given by SLHA file you have to use " /src/su/dsprep.f . 
This routine does not calculate anything because in the SLHA file al- 
ready there are all the values of the model parameters, so it only trans- 
fers the values from the file to the program. From now and on we 
concentrate only in the low energy case, that is the one in which we 
are interested. 

• Model setup 

The routine dssusy contains two subroutines: " /src/su/dsmodelsetup.f , 
and " /src/su/dsprep.f. We have already explained the content of the 
latter, that only transfer the correct values of the parameters in all the 
variables used in each subroutine. Dsmodelsetup is more important, 
because it calculates these values step by step. This is obtained using 
many subroutines: 

Constant calculation : the related subroutine is " /src/su/dssuconst.f, 
that calculates the values of the couplings, of the CKM matrix, of 
the Yukawa parameters, . . . 

Spectrum mass : the spectrum is calculated calling " /src/su/ 

dsspectrum.f that calls many subroutines: " /src/su/dssfesct.f'for 
the sfermion masses, " /src/su/dsneusct.f'for the neutralino masses, 
" /src/su/dschasct.f'for the chargino masses, " /src/su/dshigsct.f " 
for the higgs masses. These subroutines also check that the re- 
sulting spectrum is physically acceptable and stop the running of 
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DarkSUSY if not. Finally, the routine checks which particle is the 
LSP. 

Vertex : the vertices are calculated calling " /src/su/dsvertx.f'that 
contains two subroutines: " /src/su/dsvertxl.f'that calculates the 
vertices related to boson vectors, neutralinos and charginos while 
" /src/su/dsvertx3.f' calculates vertices related to sfermions. 

Higgses and sparticles widths : these widths are calculated calling 
respectively " /src/su/dshigwid.f'for the higgses and " /src/su/ 
dsspwid.f'for the sparticles. 

• Experimental bounds check 

DarkSUSY checks if a certain model satisfies the experimental bounds 
on the mass spectrum directly in the routines that calculates it, as al- 
ready said, but checks the accelerator constraints on particles interac- 
tions in another routine: " /src/ac/dsacbnd.f ". This routine compare 
the results of the program with the most recent accelerator data and, 
if a model is rejected, it stops the running of the program and prints 
a message in which it explains what experimental constraint was not 
satisfied. 

6.1.2 Relic density calculation routines 

DarkSUSY computes the DM rehc density using a specifical function, " /src 
/rn/dsrdomega.f . The complete definition is: dsrdomega(omtype, fast, xf, 
ierr, iwar, nfc). Omtype and fast are inputs that have to be given by the user 
to set the type of coannihilation that the program has to take into account 
and the numerical accuracy that has to be used to perform the calculations. 
There are many options: 

• omtype = - no coann 

• 1 - includes all relevant coannihilations (charginos, neutralinos and slep- 

tons) 

• 2 - includes only coannihilations betweeen charginos and neutralinos 

• 3 - includes only coannihilations between sfermions and the hghtest 
neutralino 

• fast = - standard accurate calculation (accuracy better than 1%) 

• 1 - fast calculation: (recommended unless extreme accuracy is needed) 
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• 2 - faster and less precise method, i.e. expands the annihilation cross 
section in x 

The output are dsrdomega, that is Qh'^, xf that is x/, ierr and iwar that are 
the identificator of errors or warnings during the calculation, and nfc that 
counts the number of function calls to the effective annihilation cross-section. 
This function performs all the numerical calculation needed to numerically 
solve the Boltzmann equation, that we have described in detail in section [H?T1 
This calculation is divided in many routines, that we are going to describe. 
Prior to call the routines, dsrdomega calculates the number of coannihilating 
particles, according to omtype. The default ratio between the mass of each 
particle and that of the LSP to be added to the coannihilating particles is 1.5, 
but can be modified by hand by the user. When the routine has ended the 
scan over the particles and determined all the coannihilating ones, there is the 
call of " /src/rd/dsrdens.f " if fast is equal to or 1, " src/rd/dsrdquad.f " 
if fast is equal to 2. 

The latter simply makes an expansion of the effective annihilation rate \5A2\ 
gaining speed but losing precision. Instead, the subroutine dsrdens imple- 
ments numerical methods to solve directly the differential equation, with the 
only approximations given by the numerical precision of the method used. It 
does not directly perform the calculation: first it sets the initial value of the 
variables that are going to be used, as the temperature and the momentum, 
then calls the subroutine. " /src/rd/dsrdtab.f 'which creates a table in which 
each value of x is associated to the related value of the effective annihilation 
rate; after that it calls the subroutine which calculates the limits of the nu- 
merical integration, " /src/rd/dsrdthlim.f "; last, it calls the subroutine that 
perform the numerical calculation:" /src/rd/dsrdeqn.f. 
This subroutine numerically solves the Boltzmann equation in the form f l5.57p 
using a trapezoidal method with adaptive stepsize and termination. Its 
input is the effective annihilation cross section, provided by the function 
" /src/rd/dsrdwintp.f. These function interpolates the discrete values given 
by dsrdtab to have all the points needed by the numerical analysis. 
It remains to analyze the set of subroutines that calculate the effective an- 
nihilation rate, used by dsrdtab. The values tabulated in this routine are 
calculated in " /src/an/dsanwx.f, that, as usual, does not calculate the ef- 
fective annihilation rate, but calls the subroutines that makes the calculation 
and organizes the results. In this case the called subroutine is the first of a set 
of nested subroutines, each one performing a certain step of the calculation. 
The complete list, in order of call, is: 

dgadap : its pattern is " /src/xcern/dgadap.f. It is called by dsanwx to 
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perform the angular integration in the 9 variable. It uses an adaptive 
gaussian method to perform the integration. 



dsandwdcosy and dsandwdcoss : their patterns are: 

" /src/an/dsandwcosy.f'and " /src/an/dsandwdcoss.f. Both are used 
as integrand function in dgadap. They differ only in the fact that 

the integration variable in dsandwdcosy is changed from 6 to y = 
— — —5 — 7^ to avoid poles in cos (6) = ±1. Both multiply the 

value of the effective annihilation cross-section by 10^5, a trick useful 
for the numerical calculations. 

dsandwdcos : its pattern is " /src/an/dsandwdcos.f . It is called by dsand- 
wcosy and dsandwdcoss. It gives the numerical value of the effective 
annihilation cross-section, summing over all coannihilation processes 
cross-sections. 

dsandwdcosij : its pattern is " /src/an/dsandwdcos.f (it is a subroutine 
written in the same file of dsandwdcos) . It is called by dsandwdcos. For 
each couple of coannihilating particles this function determines their 
types (neutralino, sfermion, chargino) and calls the proper subroutine 
to calculate their cross-section. 

dsandwcosnn : its pattern is " /src/an/dsandwdcosnn.f .It is called by 
dsandwdcosij and calculates the cross-section of neutralino-neutralino 
annihilation. 

dsandwdcoscn : its pattern is " /src/an/dsandwdcoscn.f .It is called by 
dsandwdcosij and calculates the cross-section of neutralino-chargino 
annihilation. 

dsandwdcoscc : its pattern is " /src/an/dsandwdcoscc.f'.It is called by 
dsandwdcosij and calculates the cross-section of chargino-chargino an- 
nihilation. 

dsandwdcossfsf : its pattern is " /src/an/dsandwdcossfsf.f .It is called 
by dsandwdcosij and calculates the cross-section of sfermion-sfermion 
annihilation. 

dsandwdcossfchi : its pattern is " /src/an/dsandwdcossfchi.f .It is called 
by dsandwdcosij and calculates the cross-section of sfermion-neutralino 
(or chargino) annihilation. 

anstu folder routines : the pattern is " /src/anstu". All the routines that 
calculate the amplitudes of a single physical process in a certain channel 
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are stored in this folder. They are called by the various routine listed 
above. 

6.2 Modification to DarkSUSY 

We want to modify DarkSUSY routines to have the possibility to perform 
calculations in the extended model described in chapter [31 Our scope is 
to made all the changes in such a way that permits to perform calculation 
in the MSSM if we desire it. We will start describing how we modify the 
subroutines, then we will describe some examples of main programs, using the 
main programs of the original version of DarkSUSY modified in a convenient 
way. During the description we will try to use the same order of the previous 
section, to permit easy comparisons. 

6.2.1 Model defining routines 

• Common blocks 
As we have seen, in DarkSUSY the common blocks that contain the 
variables used in the model building are defined in " /include/dssusy.h 
". This file contains the variables in four subfiles. For simplicity we 
added another file, "dsanom.h " , called together with the other four. In 
this file we have stored all the new common variables. This is the list 
of all the variables we add (we write the name used in the program), 
and their definition. 

1. Masses 

mpsiano and mOano are, respectively, the mass of the stiickelino 
and of the primeino. ampsiano and amOano are two support vari- 
ables that permit to transfer the value of the related variables 
from the input to the core of the program, mzp is the mass of the 
Z'. 

2. Couplings 

gOano is the anomalous coupling (yfo- gmix(3,3) is the 3x3 matrix 
that gives the mixing among neutral vectors, fl4.73l) . cznn(6,6), 
cann(6,6) and cpnn(6,6) are the couplings of the neutralino-neutralino- 
boson vertices (14.1061) . 

3. Anomaly related variables 

b3ano, b20ano, b21ano, b22ano, b24ano, are respectively the h^^ 
of the (I3.15P : qhuano, qqano, qlano, quano, qdano, qeano, qhdano 
are the anomalous charges defined in (13.11) . 
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4. Other 

qfnew(12),afnew(12), vfnew(12), qvanonew(12), qaanonew(12) are 
vectors that contains, respectively, the electro-magnetic quantum 
number g/, the vector and axial quantum numbers QviQa (I4.90p 
and the vector and axial quantum numbers for the anomalous 
f/(l), QviQa f!4.92p for all fermions. Dsansumfanofl,dsansfanofvvvl 
will contain, respectively, the coannihilation contributions from 
flirmj) and giiSD. 

• Initializing the program 

We slightly modify the routine dsinit: we added two neutralinos, the Z' 
and their related degrees of freedom to the particles considered by the 
program. Then we fix the values of some general constant that we have 
introduced i.e. gOano = 0.1, mzp = 1 TeV, b3ano = mzp/AgOano. 

• Setting model-dependent parameters 

In the extended anomalous model we have 5 new free parameter: Ms {mpsiano) , Mq {mOano 

The value of these parameters has to be given for each model we will 

generate. So we have modified the input system of DarkSUSY. 

To give values to the model parameters by hands we have modified the 

main program to request an input for each of these 5 new parameter. 

Furthermore, we have modified dsgive_model to tranfer these inputs to 

the subroutines. Also we have added in this routine the constraints 

(13. ip and the definitions f l3.15p . 

To define the model parameters randomly we have modified the sub- 
routine random_model(), adding the 5 new parameters to those that 
have to be randomized. Because this subroutine is written in the same 
file of the related main program, we have introduced directly in the 
main program the constraints (13. ip and the definitions (13.150 . 
We are not interested in using SLHA files to introduce the model pa- 
rameters, so we do not have modified that part of the program. 

• Defining all model-dependent quantities from the inputs 

Being our model an effective model, we skip on the routines that per- 
form RGE evolution and on the routines that use SLHA files to define 
the model dependent quantities. Thus we concentrate only on dssusy. 

• Model setup 

We have modified the subroutines called by dsmodelsetup in order to 
describe the extended model: 

Spectrum mass : first of all we modify the command that checks 
if the LSP is a neutralino in such a way that it considers all six 
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neutralinos and not only the four of the MSSM, then we add the 
definition of axinofrac as the fraction of anomalous neutralinos 
in the LSP. We have also modified the subroutines dsneusct and 
dshigsct, that calculate the masses of neutralinos and higgses. In 
dsneusct we have rewritten the mass matrix, substituting fl2.6p 
with fl3.45p and changing all the cycles lenghts from 4 to 6 to take 
into account the two extra neutralinos. 

Modifying the higgs masses calculation is not straightforward. 
The reason is that the latest version of DarkSUSY does not cal- 
culate these masses with inner routines but using the external 
program feynhiggs. We want to modify tree-level values for higgs 
masses, implementing f l3.3ip . To achieve this goal, we have to 
modify a feynhiggs subroutine. DarkSUSY stores the external 
programs that it utilizes in the folder " / contrib " . In this folder 
there are many version of feynhiggs. The one used by DarkSUSY 
is 2.6.3. The subroutine in which there is the tree- level higgses 
mass calculation is " /contrib/FeynHiggs-2.6.3/src/Main/Para.F 
" . Here we have changed the formulas according to (13.311) . 

Vertex : we had to change the vertices according to the various equa- 
tions of section I4l3l The main routine that implements the vertices 
in DarkSUSY, as we have already said, is dsvertx, that contains 
dsvertxl and dsvertx3. These are the files that we have to mod- 
ificate. 

We start from dsvertxl. First we introduce the definitions of the 
variables related to the new vertices: so we give the definition of 
gmix according to (14. 73 p . the definition of vfnew, afnew, qvanonew 
and qaanonew according respectively to the SM definitions [32j for 
the first two, and to (I4.92p for the others two. After that it starts 
the vertices definition part of the subroutine. First of all there are 
the MSSM vertices, already defined in DarkSUSY but that have 
to be modified to describe our model. So we modify the higgs- 
higgs-vector boson vertices according to (14.1 16p and the fermion- 
fermion- vector boson vertices according to (14.930 and (14.970 . Then 
we extend the length of the cycles in which there is the definition 
of the vertices in which appear neutralinos from 4 to 6, to take 
count of the two extra neutralinos of our model. After that we 
added the new vertices: 7-neutralino-neutralino, Zo-neutralino- 
neutralino, Z'-neutralino-neutralino according to (I4.104p . 
In dsvertx3 we have modified the sfermion-sfermion-vector boson 
vertices according to (I4.93P and (I4.97P and the sfermion-sfermion- 
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higgs vertices according to (I4.9UI) . 



6.2.2 Relic density calculation routines 

We have not changed the subroutines that perform the numerical integration 
of the Boltzmann equation, as the equation obviously does not change intro- 
ducing the anomalous U{1). As we have seen, the changes are the presence 
of new interactions and of new particles, as well as some modifications to the 
MSSM interactions. DarkSUSY already contains all MSSM interactions, so 
their modifications, listed in the previous section, are already taken in ac- 
count, while we have to implement the new interactions f l4.122p and fl4.125p . 
We want to keep separated the routines that we have written to add the 
contributions to the cross section that were not present in the MSSM and 
consequently in the DarkSUSY package. So we have created and added to 
the makefile the folder " /src/anano ", where we have stored the new sub- 
routines that we are going to describe. 

• dsanscalarproduct 

This subroutine contains the definitions of the kinematical variables 
fl5.79p . It calls dsankinvarl, a modified version of dsankinvar, a Dark- 
SUSY subroutine that contains the definitions of the kinematical vari- 
ables already used by the program. 

• dsankinvarl 

It is identical to dsankinvar except for the definition of the mass of 
the exchanged particles, that we skipped because we have not common 
vertices for all the vector boson, as in the original program, so we have 
not the need to distinguish among them. 

• dsansumfanof 

This subroutine resets the value of dsankinvarl 0, and it sums the 
contribution of (14.1220 over all fermions of the SM. 

• dsansfanofvff 

This subroutine calculates ( I4.122p for every couple of neutralinos and 
SM fermions. 

• dsansfanofvvv 

This subroutine calculates fl4.125p for every couple of neutralinos. 

^It is necessary because the kinematics changes when we change the studied process 
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The only change we made in the original DarkSUSY relic densiy routines are 
in dsandwdcosnn, where we add the call of dsansumfanof and dsansfanofvvv 
and a final sum of a 10^^'' that is a rough estimate of the contribution of 
the annihilation of an LSP mostly composed by stiickelino. This is needed 
to avoid the possibility Weff = that causes the failure of the numerical 
algorithm to solve the Boltzmann equation. Because this contribution is an 
estimate, it makes the result of a calculation in which there is only the anni- 
hilation of a stiickelino-type LSP absolutely imprecise, but as we have seen 
in section 14.11 this contribution cannot ever explain WMAP data, while if 
we have coannihilations this contribution is negligible, so this approximation 
does not modify the physical results. 

6.2.3 Main programs 

We have written two main programs. The first (anomssmS) has the scope to 
study random models, while the second (dsmain) permits to study a specific 
model. 

Because anomssmS has the purpose to study a complete randomly generated 
model we have modified the DarkSUSY subroutine random_model adding 
the five parameters from the anomalous extensions to be generated along the 
seven parameters from MSSM already implemented. The new parameters 
are Qq, Ql, Qhu^ and Mq, implemented respectively as qqano, qlano, 
qhuano, mpsiano, mOano. We add an user interface to decide if we want to 
perform simulation in the MSSM or in the miaUSSM, if we want a bino- 
type or wino-type MSSM LSSP, to introduce manually the range in which 
each parameter have to be generated and the number of models that will be 
generated. Then we add a command to write all the twelve parameters in 
a file that we had already open in the main program. In this way we have 
stored the input parameters of each model. 

Because the models will be randomly generated, we have not control on the 
mass gap between LSP and NLSP. It will be a result that depends on many of 
the model parameters. However we want to have the possibility to study the 
results with a "controlled" mass gap. So we add checks after the relic density 
calculations that permits to store the results in different files if the mass gap 
is in the following intervals: - 5%,5 - 10%, 10 - 15%, 15 - 20%, > 20% plus 
a file in which are stored all the results, that can be useful if we want an 
ensemble vision. 

The calculation part of this main program is very simple: first we randomly 
generate the parameters of the model with the subroutine random_model, 
then we build the model with the subroutines dsmodelsetup and dssusy, then 
we check the accelerator bounds and the constraints on the neutral mixing 
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(subroutine acbounds), then we perform the rehc density calculation (with 
dsrdomega) only if the LSP satisfies certain conditions (for example we can 
request that it is a combination of stiickelino and primeino). After that the 
relic density results are stored (with other important quantities calculated, as 
Xf, the spectrum,...) in the file related to the mass gap. We store separately 
the models that satisfy WMAP from those that doesn't. 
The second program, a modified version of dsmain, a program already present 
in the current DarkSUSY distribution, permits to the user to enter by hand 
the values of the twelve parameters that define the model. This is made 
using explicit write and read commands for each quantity. After that we call 
the already modified subroutines dsgive_model and dssusy, then we call the 
same subroutines of the previous case, acbounds and dsrdomega. Because 
this program deals with one user-defined model it is not necessary to open 
files for many possible mass gaps, because we will have only one result. 
Note that, if in anomssmS we choose the lower and the higher extremes of 
the parameters range to be equal and if we set the number of models to 1, 
the program works in the same way of dsmain. We can use this property 
to perform simulations in which we keep constant some parameters to study 
the dependence of the model from the others. 

6.3 Simulations and results 

As for the study of the no mixing case, we perform our simulations for the 
different situations of a bino-higgsino NLSP and a wino-higgsino NLSP. In 
both cases the LSP will be a random mix of stiickelino and primeino. 
We have started performing an ample scan (number of model > 10^) in a 
wide area of parameter space, without any constraint. We have utilized these 
results to find the regions of parameter space in which there is a model with 
relic density that satisfies or is near to the WMAP results. 
Tipically, in each region we have found, the mass gap between LSP and NLSP 
of the models that satisfy WMAP data is almost constant, but from region 
to region it can change from 1 2% up to 25%. 

To perform a more accurate study, we have chosen sample region to be ex- 
plored more carefully for each of this mass gap ranges: -v- 5%, 5 ^ 10%, 
10 ^ 15%, 15 ^ 20%, > 20. Then we have used our modified DarkSUSY to 
scan these regions keeping all but two parameters constant. 
In this way we have obtained the dependence of the relic density on each 
couple of parameters in regions in which WMAP data are satisfied. Because 
now we are interested in well defined regions of parameters space we have 
used less (~ 10^ instead of ~ 10^) number of models. 
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Figure 6.1: Plot of relic density versus Qh^ and Qq (left, bino-liiggsino 
NLSP) and versus Qh^ and Ql (right, wino-higgsino NLSP) 

6.3.1 General results 

In this section we want to list some results that are valid for both types of 
NLSP. 

First of all, given the constraints on the neutral mixing described in [3S] , we 
want to found their effects on the model parameters. The results are showed 
in figure [6?H for two sample models. These plots show the general result: 

Constraints on the mixing ^ — 1 ^ Qhu ^ 1 (6-1) 

As we can see from the figure 16. 1[ the contribution of the anomalous charges 
to the relic density is extremely small. 

6.3.2 Bino-higgsino NLSP 

If the NLSP is mostly a bino-higgsino we have a two particle coannihilation. 
We want to study the relic density with respect to the parameters that mostly 
contribute to the LSP mass {Ms and Mq), with respect to the parameters 
that mostly contribute to the NLSP mass (M2 and /i) and also with respect 
to the other MSSM7 parameters. The anomalous charges Ql, Qq, as 

we have seen in the previous section, does not contribute in a significant way 
to the relic density. 

We studied the relic density with respect to Ms and Mq for different values 
of the mass gap. The results are shown in fig 16. 2[ As we can see, changing 
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Figure 6.2: Plot of the relic density of the LSP vs the stiickelino and primeino 
masses in the case of bino NLSP. The first plot shows the case of mass gap 
5 %. The second plot is a zoom of the first in the region between 20-140 
GeV for stiickelino mass and 500-850 GeV for primeino mass. The third plot 
shows the case of mass gap 10 % in the same region of the second plot. All 
the masses are expressed in GeV 

the mass gap from 5% to 10% greatly reduces the region in which WMAP 
[37] data are satisfied. Only in the region showed in the second image the 
relic density satisfies WMAP data for both cases. 

We have also studied the relic density with respect to M2 and /i. The results 
are showed in the figure 16. 3[ The image is roughly symmetric with respect 
to the line /i = 0, also showing that near this line we do not satisfy WMAP 
data, while in the TeV range for fi we satisfy WMAP data in each point that 
is not forbidden by experimental constraints. Furthermore, we can see that 
the region of low M2 is physically unacceptable. 

The figure obviously refers to a sample model, but the fact that the relic 
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Figure 6.3: Plot of the relic density vs M2 and /x, MSSM parameters, for 
mass gap 10% in the case of bino NLSP. The masses are expressed in GeV 

density behaviour is symmetric with respect to /i = and that for low M2 
there are violation of experimental constraints is a general result. 

6.3.3 Wino-higgsino NLSP 

If the NLSP is mostly a wino-higgsino we have a three particle coannihila- 
tion, because the wino is almost degenerate in mass with the lowest chargino. 
Again in this case the anomalous charges Qhui Ql, Qq do not contribute 
significantly to the the rehc density. So we performed the same study de- 
scribed in the previous section. 

The results for the dependence of the relic density on M5 and Mq are very 
similar to those of the previous section, so we will only show one of them 
as an example in fig. 16. 4[ Furthermore can be very interesting to show the 
presence of a phenomenon called funnel region. It is not new: it is already 
predicted in the MSSM. Substantially it consists in a sort of resonance in the 
case in which the relation 2Mlsp ~ holds. So there is a peak of the 
relic density versus the mass of parity odd higgs boson Aq. In the MiAUSSM 
we have a different LSP, the stiickelino-primeino mix, but the NLSP is the 
old MSSM LSP, whose interactions contribute greatly to the relic density so 
we expect the funnel region phenomenon to appear in the same way of the 
MSSM. As it is showed in figure 16.51 the funnel region appears as expected. 

We can also show for completeness our result for the dependence of the 
relic density on and M2 (fig. 16.61) . 
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Figure 6.4: Plot of the relic density vs Ms and Mq for mass gap 10% in the 
case of wino NLSP. The masses are expressed in GeV 
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Figure 6.5: Plot of the relic density vs M^o and mo, MSSM parameters, for 
mass gap 10% in the case of wino NLSP. The masses are expressed in GeV 




Figure 6.6: Plot of the rehc density vs // and M2, MSSM parameters, for 
mass gap 10% in the case of wino NLSP. The masses are expressed in GeV 
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Chapter 7 
Asymmetry 



7.1 Introduction 
7.1.1 PDFs 

The study of the asymmetry at the LHC imphes deahng with hadronic pro- 
cesses. To calculate cross sections that involve hadrons we need the Partonic 
Distribution Functions (PDFs) [IQ]. These functions give the distribution 
of momenta of the elementary constituents of a given hadron in the parton 
model. In this model the elementary constituents of the hadrons, called par- 
tons, are point-like and almost free particles. Because QCD is an asymptotic 
free theory [40j this assumption is accurate for high momentum transfer (that 
is the so called Deep Inelastic Scattering Region). So this model holds for 
momentum exchanges much larger than the mass of the particles. We are 
surely in this situation at the LHC so we will use this model in the follow- 
ing sections. Parton is a generic term that indicates the constituents of the 
hadrons. So a parton can be a quark or a gluon. 
To introduce the PDFs let's consider a generic process: 

A + B ^c + d + X (7.1) 

where A and B are the initial hadrons, c and d are the final particles in which 
we are interested and X represents the hadronic remnants. The elementary 
subprocess associated to the hadronic process is: 

a + b^c + d (7.2) 

where a and b are elementary constituents of the hadrons. This cross section 
can be calculated with the usual Feynman rules derived from the Lagrangian. 
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We need the PDFs to obtain the cross section (17. ip from the elementary one. 
Infact (17. ip is given by the convolution of fl7.2p with the PDFs: 

da{AB cdX) = (7.3) 




where: 

• a{AB —7- cdX) is the cross section of the hadronic process 

• a{ab — )■ cd) is the cross section of the elementary process 

• /(a/b)/(A,B) (a^(i,2), Q^) is the PDF that gives the probability to have the 
parton a/b in the hadron A/B with momentum fraction Xi/2 and with 
momentum transfer (this is sometimes omitted for simplicity). 

Being Xj a fraction of the total momentum of the hadron that contains the 
related parton, we have necessarily: < a;j < 1. 

The PDFs are obtained using experimental data, because the present numer- 
ical calculations in QCD are not able to give satisfactory results. So there 
exist many sets of PDF, given by different research groups worldwide. For 
the calculations that we are going to perform in the next sections we will use 
the PDFs defined in [H] and 02]. 



7.1.2 Asymmetry definition 



The asymmetry |44j is a quantity that, given a certain particle interaction, 
measures the difference between the number of final particles emitted in one 
direction ("forward") and in the opposite direction ("backward"). The direc- 
tion is usually chosen with respect to that of the initial particles: we define 
as forward the semisphere identified by the beam direction, while backward 
is identified by the opposite direction. The basic equation is 

F-B 

Afb 



where 



Explicitly: 



A 



FB 



F + B 

F = ^ of particles going forward 
B = ^ of particles going backward 

f dXldX2 J2g fq{Xl)fq{X2){F - B) 
Jdx,dX2Ej<li^l)fqi^2){F + B) 



(7.4) 

(7.5) 
(7.6) 

(7.7) 
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Figure 7.1: The SM LO ApB prediction as a function of the dielectron in- 
variant mass for uu — )■ e+e^, dd — )■ e~^e~ , and pp — >■ e+e" 



with: 



^ , „da(cos6, s) „ , ^da(cos9, s) , 

dcose—r — B= / dcose—^ — -f^ 7.8 
d cos dds j_i d cos 6ds 

In these formulas xi and X2 are the standard partonic variables, so the integral 
over them goes from to 1. fq{xi) are the partonic distribution functions 
(PDFs) of the quark q with respect to the variable Xj. cos 6' is the angle of 
emission of the electron of the final state in the CM frame, ^"^^^"^^f^ jg the 

a cos yds 

differential cross section in the CM frame. 

The asymmetry has been calculated and measured for the SM with good 



precision. A sample of the results obtained |13] are in figure 17.11 It has been 
used for many phenomenological purposes, such as giving constraints on the 
mixing angles of SM, measuring the couplings, setting constraints on the new 
physics beyond the SM. 

At the LHC it can be used in the same ways as well to discover new physics. 
In the following we aim to calculate the asymmetry at the LHC to impose 
constraints on the free charges of the MiAUSSM and to give a method that, 
if an asymmetry different from the MSSM predicted one will be measured, 
permits to check the consistency of our and other models that extend the 
MSSM. 
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7.2 Four Asymmetries at the LHC 



In e^e^ or pp colliders, in which the initial state is asymmetric, the asymme- 
try is naturally different from 0. However at the LHC the total asymmetry 
is zero since the initial state, that is proton-proton, is symmetric. The el- 
ementary subprocess qq — )■ // is instead asymmetric. This property is lost 
when we integrate over the full space of partonic and kinematic variables to 
calculate the asymmetry of the pp process. However if we do not consider 
the total parameter space performing suitable cuts we can obtain a non-zero 
value for the asymmetry. Consequently we can have different definitions of 
asymmetry at the LHC changing the choice of the cut. In the following we 
present four definitions of asymmetry viable at LHC, each one associated to 
a different cut of the parameter space. These can be presented in couples: 
the Forward-Backward asymmetry with the One-Side asymmetry and the 
Central asymmetry with the Edge asymmetry. 

The definition of the "Forward-Backward" asymmetry is [52], j45] . |46] . 

mm- 





> 


Yf\)-cri\Yf 


< 


Yf\) 


\Y\>\Y^ut 




> 




< 


Yf-l) 



(7.9) 



where a is the cross section after having performed all the integrals and Yy, Yj 
are the pseudorapidies of the outcoming fermion or antifermion, respectively. 
In the Centre of Mass (CM) frame their expressions are: 

>///-=- log (tan (7.10) 

As already said in the previous section 6f/f is the angle of emission of the 
electron/positron of the final state in the laboratory frame. To perform the 
calculation we need to write 17.91 with respect to 6 that is the angle of the 
fermion with respect to the beam axis and is the simplest angular variable 
that we can use. This angles are represented in figures UT^l and 17^51 

So we have to rewrite the conditions \Yf \ ^ \Yj\ with respect to 6. In the 
CM frame 6f = 6f + iT, so we have: 




(7.11) 
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Figure 7.2: Representation of the process qq — )■ // in the CM frame, with 
exphcit indication of the angle that we have defined as 6 



Lab frame 




Figure 7.3: Effect of the boost from the CM frame to the lab frame. Note 
that the final particles is squeezed in the beam direction while the final 
antiparticles is squeezed in the transverse direction 
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Rewriting the related pseudorapidity we obtain: 

'6-. 




= — log ycot 
-Yf (7.12) 



Our conditions \Yf \ ^ \Yj\ become: 

tan (^1^ > cot (t}j^_l<0^<l^ cos9f > (7.13) 
6f\ ( 6f\ TT „ Stt 



tan (^^J < cot (^^J ^ - < 0^ < — ^ cosOf < (7.14) 

So we can identify the forward direction as the region with cosOf > and 
the backward direction as the region with cosOf < 0. From now on, we will 
use only the angle 6f for the calculations in the CM frame, so we will call it 
simply 6. Later we will use the expressions 6 j/j to indicate the angles of the 
outcoming fermions in the lab frame. 

Having just identified the forward and backward regions, we use the definition 
in the CM frame (TO 

dcos9—^ -f^ B= / dcosO—^ -f^ (7.15) 

g d COS Ods J_i rf cos yds 

These expressions identify the forward and backward part of the cross section 
associated to the fundamental subprocess qq — )■ e~^e~ . The relation (17. 9p can 
then be rewritten 

^^^jx,dx,Y.,U^i)h{^2){F + B) ' 

where Ccut is the integration domain and it is dependent from the cut > 
\Ycut\- As we have already defined fg, fq are the parton distribution functions 
(PDFs) associated to the quark q or to the antiquark g; xi and X2 are the 
fraction of the proton momentum carried by the quark or the anti-quark, 
respectively. Xi and X2 are not the best variables to perform the integration, 
because the standard calculation of Feynman amplitudes gives us the differ- 
ential cross section in terms of kinematic variables, such as s and cos 9. So 
we want to write the integrals in terms of s and the rapidity Y, with the 
definitions: 

s = Sx^X2 Y = \\og('^) (7.17) 
2 Vx2/ 
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S = (14 TeVf at the LHC. This relation for Y is vahd only in the CM 
system and in the ultra-relativistic limit. In fact the general definition for 
the rapidity is 



E + p, 



- log ( — — 

2 ^ VI - 



1 + Vrel j 



Vrel 



(7.18) 



with Vrel = Pz/E being the relative velocity between the two interacting 
particles. 

Using the new variables the formula for the asymmetry can then be written 
in the following form 



A 



RFB 



J dsJ 


r+oo 

^ ^cut 


dY+ f ^^^'dY 




[s,Y))f,{x2 


>,Y)){F- 


B) 


J dsJ 


r+oo 

^ ^cut 


dY + r"^'="* dY 

J —oo 


J2q fqi^l 


:s,Y))fgix2 


:s,Y)){F + 


B) 



(7.19) 

J is the Jacobian associated to the change of variables. Using 17.171 we can 
obtain 



Xi 



Y 



So the Jacobian is 
J = 



dgXi dyXi 
dsX2 dYX2 



X2 



1 



2VsS 



-Y 



2VsS 



-Y 



1 

S 



(7.20) 



(7.21) 



As we can see, in formula 17.191 the cut used in Ajips is explicit and consists 
in excluding the region —Ycut < Y < Y^ut- We will see in the following that 
this domain of integration can be further reduced limiting the range of the 



variable s. Now we focus on the integral dY: if we perform the change 

of variable Y — )■ — F, because of the relation in fl7.17p . this corresponds to the 
exchange Xi ^ X2 and consequently to the exchange of what is forward with 
what is backward {F -H- B). By summarizing, the final effect of Y —Y is 



dYj2f,i^i)M^2){E±B) 



dYj2f,ix2)f,ix^)i±F + B) 
q 

(7.22) 



Finally the FB asymmetry in fl7.9p can be rewritten in the useful form 

/ dsJ dY {s, Y),X2{s, Y)){F- B) 



A 



RFB 



with 



/ dsJ dy E, fUMs, Y),X2{s, Y)){F + B) 



fqq{.Xl,X2) = fq{Xl)fq{X2) ± fq{Xl)fq{X2) 



(7.23) 



(7.24) 
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Another type of asymmetry defined in [IH] , [SO] is the " One-side" 



A 



o 
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Yf\)-cr{\Yj 


< 




\Pz,ff\>\Pz,cut 




> 




< 


yf\) 



t\ 



(7.25) 



the only difference with the previous one is in the cut, no longer performed 
on the rapidity but on the longitudinal momentum along the beam axis, 
Pz- We can apply the same steps described before with now J_^'''''* dpz 

instead of /_^™' dY (the exchange is now — t- —pz for xi ^ X2). Because 

Pz = ^{xi — X2) the change of variable is now 



Xi 



(pz + VpfTs) X2 = ^(^-Pz + ^pI + s 



(7.26) 



with obviously a different definition for J: 



J 



dsXi dp^xi 

dsX2 dp^X2 



2^5(s+pi) 
1 



J-(l J 

J-f-l -I 2^ 



This leads to the result: 

J ds l^^2t^P^'^T.qfqgi^lis,Pz),X2{s,Pz)){F - B) 



A 



o 



JdsJ °^ dpzJj2qfqqi^li^^Pz),Ms,Pz))iF + B) 



S\^s+pI 
(7.27) 

(7.28) 



The final two definitions of asymmetry, which explore complementary regions 
in the space of angles, are the "Central" [51], [52], [53], [M], [SS] and "Edge" 
[56] asymmetries 



Ac{Yc) 
Ae{Yc) 





<Yc)-aj{\Yj 


<Yc) 




<Yc) + aj{\Yj 
>Yc)-aj{\Yf 


<yc) 

>Yc) 




>Yc) + (yj{\Yj 


>yc) 



(7.29) 
(7.30) 



In these two cases the cut is not performed on the rapidity nor on the lon- 
gitudinal momentum but on the angle (therefore on the limits of integration 
for F and B). These asymmetries are defined in the lab frame, because the 
Lorentz transformation from the CM frame "squezees "the final particles. So 
in the laboratory (Lab) frame the outgoing e~ and e"*" have no longer the 
same direction: the angles ±^ with respect to the z axis for the respec- 
tively in the CM frame are replaced by the different 6^ and 6^^ . So the 
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Figure 7.4: Graphical representation of the Central and Edge regions in which 
we have to integrate to obtain the Central and Edge asymmetry, respectively 

explicit form of the definitions fl7.29p in the lab frame is: 

From this expression we can understand the different names of the two defi- 
nitions: the "Central" asymmetry has the angular integration on the central 
angular region, the "Edge" in the complementary region. This is well shown 
in figure [731 As we have already seen the boost given by the Lorentz transfor- 
mation from the CM frame to the Lab frame "squeezes" the e~ mainly in the 
forward and backward directions, while the e"*" are left more abundant cen- 
trally [51J: according to the previous definitions Be > Fq and consequently 
we expect negative values for the Central asymmetry (in all the other three 
cases we have always positive values). 
Their expressions in terms of F and B are: 



^'^ f dXldX2 Y.q fq{Xl)fq{X2){Fc/E + Bc/e) 



Because we have calculated the cross section of the processes in the CM frame 
we have to calculate how the definitions (17.311) appear in the CM frame. In 
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this frame the differential cross section for e"*", here denoted with the tilde, 
is obtained from that of the e~ by simply inverting the sign of the cosine of 
the angle 

rlrr rin 

(7.33) 



cos S— >— COS ^ 



d COS 9ds d cos 9ds 
So the previous definitions f l7.3ip for F and B in the Lab frame assume the 
following form in the CM frame: 

Fc = f^f -^^d cos 9 Be = fT''' TT^d cos 9 

J— 1 dcosOds J — 1 d cos Ods 

Fe = ftr^ :r^^dcos9 Be = ft^ -r^^d cos 9 (7.34) 

^ J J {cut) d COS Ods ^ J J {cut) d cos Ods ^ ' 

where /(x) is simply the relativistic function which relates the cosines of 
angles in the Lab frame to those in the CM frame, i.e. cos^*"^ = f{cos9^°'''): 

f{x) = 2 , . ^ , 7 = 2 ; Vrel = ] 7.35 

-7^ + (-1 +7^)x^ ^1 - v^^i X1 + X2 

The calculation of this function is a simple application of special relativity. 
We start from the equation that gives the relativistic velocity in the lab 
frame in function of the relative velocity with respect to the CM frame and 
the angle between the initial and the final fermion in the CM frame (this is 
the angle that we have used for the first two definition of the asymmetry at 
the LHC, 9) p]: 



(„, \ f sinf^ 

l^/M' 7(1+1)2 cos6») 
/ \ v{l+cosO) 

y^f)--^ ~ (1+1,2 COS e) 



(7.36) 



In this formula z is the beam direction, x is the orthogonal direction, u is the 
velocity of the electron in the lab frame and v is the relative velocity between 
the lab and the CM frames. Now we want to calculate the relation between 
cos9f and cos 9. 

tan9f = = -j^^ = A(cos9) 

J {uf)^ 7(l+cose) V ^ 2n a2 2n 

■ n ^^ ^ 1 — COS 9t = A cos 9t 7.37 

, r, sm0f -v/l-cos^ef / / V / 

tanfe*/- = — f- = a — 

J cos f cos f 

Solving the equation on the right we obtain: 

cos9f = ^ = (7.38) 

1 _|_ 1— cos2 
^ 72(l+cose)2 

If now we invert this relation we have: 

cose^ "r'^l^yt^, (7.39) 

-72 + (-1 + 72)cos2 6'/ ^ ^ 
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that is the formula (17.351) that we wanted to demonstrate. 

The definition of asymmetry for the Central and Edge cases can therefore be 

presented as follows: 



with 



J 



dsXi dx^xi 

dsX2 dx^X2 



Sxi 



Fc/E — Bc/E 




1 

25x1 



Sxl 



2xiS 



(7.40) 



(7.41) 



Instead of the angle cut it is convention to write the cut performed on these 
two definitions of asymmetry using the related pseudorapidity 



Yr 



c 



log tan(r"72) 



(7.42) 



where = arccos(ctit) This variable (usually denoted with rj in literature), 
must not be confused with the cut on the rapidity Y previously used for the 
RFB definition of asymmetry. 



7.3 Cross-section 

As we said earlier we are interested in the cross-section in the CM for the 
process pp — )■ e+e~. This is the convolution of the elementary cross-section of 
the process qq — )■ e^e" with the PDFs of the proton. So we have to calculate 
the cross-section of qq — ?■ e^e~ for a general Drell-Yan interaction. This 
means that we can have the product of diagrams in which can be exchanged 
the 7, Zq and Z' , showed in fig. 17.31 Thus we can have 6 possible terms in 




q 



Figure 7.5: Feynman diagrams for the Drell-Yan process in the MiAUSSM 
the amplitude: 77, 'jZq, jZ', ZqZq, ZqZ' and Z'Z'. We will sum them all: 

\Mnq) = l\ E alglM^b (7.43) 

a,b=-y,Z,Z' 
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where the fractions | and | come out from the averages over spin and color, 
Qq is the couphng associated to the Z'. We fix its value to 0.1. Mab is the 
amplitude of each process divided by the couplings. 

To calculate Mab we have to write down the Feynman amplitudes divided by 
the couplings for the processes in figure 17.31 for the exchange of any of the 
three possible vector bosons: 

= z [uf^,{Cl, - l,Cf^a)vf] — ^ [vMCla ' l^C^J^,] (7-44) 

^ a 

^J/q a' ^J/q a quautum uumbers of the quark q or the fermion / with 

respect to the interation a, that goes from to 2 and represent the U{iy, 
[/(I) and SU{2) respectively. Their values are listed (with the electron as 
the fermion /) in table 17.11 From now on we use the subfix e instead of / 
because we are focusing on electrons in the final state, but the formulae can 
be used for other fermions simply using the appropriate quantum numbers. 
Its square modulus gives: 

Ma, = (s-M^){s~M^) ^''^^^^ + ^eh,{Cl, - l,Ct,a){h ' ^/) (7-45) 

Calculating the products, using the properties 757^ = —7^75 and 7I = 1 and 
the fact that irtg = rrig, rng = rrig the first trace becomes: 

tr[{^, + m,)7,(C]:, - 75C^^J(^/ - mj){Cl, + 75^^^)7.1 = 

ClaCibtr[{^,l,Pel5lu] + (7.46) 
CtCl,btr[^M-l^)Pe7>^] + 

In the same way for the second trace we have: 

ir\{i, - "^q)l^(Cla - l,Cla){i, + m,){Cl, - Ci^Y\ 
ClaClbH^.l^Pql" - m,m,-7'^7l + 

ClaCtifrWqrPql^Y] + (7.47) 

CtaClbtr[i>,l''{-l,)PqY] + 

CtaCtbtr\^{l'PqY + m,m,7'^7l 



95 





f Z' 




f = u,c,t 


Qq + QhJ2 


-QhJ2 


f = d,s,b 


Qq - QhJ2 


QhJ2 


f = e,iJ,,T 


Ql-QhJ2 


QhJ2 



Table 7.1: Vector and axial quantum numbers of the SM fermions with re- 
spect to the Z' 



Multiplying the two contributions and using the formulae fl4.20p . fl4.2ip . 
(1122D, (i2SD and f lQIj) we obtain the final expression: 



M, 



ah 





{s - ml) {s - ml) + (r^ ma Tf, mf,) 




{s - mlY + {TamaY 





X 



2 ml m^ {C^a C^e,b ~ ^'e,a Cl,b) {^t,a ^tb " ^g^a^'Tf') ~ "^e {Pq ' Pq) 



X 



X 



ipta Ce,b ~ Ce,a Ce,b) (^^a Cq,b + ^q,a ^^fe) + [Pq ' Pe) [Pq ' Pe) 
{pl, Cla + Cla Cl,) [Cf,, Cla + C^, + [p, ■ Pe") [p, " Pe) X (7.4^ 

^e,a + ^^,a ^I,b) {^q,a. ^q,b + ^q,a ^q,b) ~ {Pq ' Pe) {Pq ' Pe) X 

i^e,b ^e,a + ^e',b) {^tb ^q^,a + ^q,a C'^fe) + {Pq ' Pe) {Pq " Pe) X 

f/^A ^A , \ I r<A (^A , ' 

\^e,a ^e,b + ^ e,a ^ e,b) q,a ^ q,b + ^ q,a ^q,b. 

( r~iA r~iA , r~iV r~iV \ ( r^A r~iA r~iV r~iV \ 
\ye,a ^e,b "r ^e,a ^e,b) \^q,a ^q,b ^q,a ^q,b) 



ml {pe ■ Pe) 



In this expression we have substituted the divergent expression of the propa- 
gator, — with the regulated one where is the width of the 

boson a. Nab is a multiplicity factor that is equal to | if the exchanged vector 
bosons are identical and is equal to 1 if they are different. The Cs are simply 
the vector and axial quantum numbers related to the vector bosons: for the 
7 and the Zq they are the usual SM quantum numbers that can be found 
in [32], while the vector and axial couplings of the extra U{1) have been 
previously calculated in fl4.92p and are showed in table 17.11 The differential 
cross section can be found multiplying for the usual kinematic prefactor and 
summing this result over the contribution of the 6 possible initial quarks: 



9V 



dsd cos 9 



CM 



E 



Pe 



32 IT SPq 



m'iq) 



(7.49) 
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To obtain its explicit expression we have to write the scalar products with 
respect to s and 6. Using the formulae fl5.79l) we obtain (in the CM frame): 



(PgPg) 






(PePe) 


= E,E,+ \p,\^ 




(PqPe) 


= EqEe - \Pq\\Pe\ 


COS 9 


(PqPe) 


= EgEg + \Pg\\Pe\ 


cos 9 


(PqPe) 


= EqEe + \Pq\\Pe\ 


cos 9 


(PqPe) 


= EqEe - \Pq\\Pe\ 


cos 9 



(7.50) 



Because the couples of incoming and outcoming particles have the same 
masses we can use the definition of s to calculate all the quantities in this 
expression. In the CM frame we have: 

ml = {El - = {El - \pe\^) ^Ee = E, (7.51) 

so: 

S = {Pe+ Pef = {Ee + ^e, 0)^ = {2Eef ^ Ee = ^ (7.52) 

\Pe\ = C^-SS) 

Substituting e with q we obtain the corresponding relations for the quarks. 



7.4 Asymmetry calculation 
7.4.1 Optimized asymmetry 

As we have explained, each definition of the asymmetry at the LHC contains 
a certain cut on the parameter space. The cut acts as a variable of the 
asymmetry so it has to be fixed in order to obtain a numerical value from 
the calculation. The problem is to choose the optimal cut. One can be 
tempted to use the value for which the asymmetry is maximum, but this is 
not the best choice. 

To find the best cut, we have to define the significance: 

Sig = A^C Gtot (7.54) 

with £ = 100 fh~^ as the integrated luminosity at LHC and a tot the total 
cross section of the process, sum of the forward and backward contributions. 
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So the optimal cut is the value of Y^ut, Pz,cut or Yc for which the significance 
is maximum. 

It depends mainly on the properties of the PDFs and on the mass of the 
Z'. In this work we set Mz' = 1.5 TeV. Because the best cut is nearly 
independent on the properties of the Z' coupling to fermions [57], we can 
find it fixing the free charges Qi associated to the extra U {!)'. Obviously the 
values of the asymmetry and the significance are dependent on these charges: 
we will concentrate on them in the next part of our analysis. 
Furthermore it is useful to restrict the range of s of the studied region: the 
standard choice is to use the "on peak" and the "off peak" regions, defined in 
the sequent way: 

"on peak" ^ Mz' - 3Tz' < M^+e- < Mz' + 3Tz' (7.55) 
"off peak" ^ '^Mz' < M^+e- < Mz' - 3Tz' (7.56) 

We will use only the "on peak" region for simplicity, but our calculation can 
be perfomed in the same way in the "off peak" region. 

As said, we have performed calculations with fixed charges in order to find 
the best cuts that then we will use to study the dependence of the asymmetry 
from the charges. Our choice has been Ql = 1, Q^^ = 1/2 and Qq = 3/4. 
The results for the best cuts, that are obtained from the figures in the fol- 
lowing, are: 

As a check, because the cuts have to be nearly independent from the charges 





Arfb 


Ao 


Ac 


Ae 


Best cut 


Yff = 0.4 


p'^^j = 580 GeV 


Yc = 0.8 


Yc = 1.4 



Table 7.2: Best cuts for the on-peak e~^e asymmetries 



[57] , we have calculated the cuts after changing the value of Qhu^ Qq and Ql, 
obtaining the same results just listed. 

FB asymmetry and significance with fixed charges: the search of 
the peak for the significance 

The figure 17.61 shows the dependence of the on peak asymmetry and the sig- 
nificance versus the cut Ycut- It is clear that the significance has a maximum 
around Ycut = 0.4. Using this result we have obtained the figure ITTTt in wich 
we keep fixed the cut and study the dependence of the asymmetry and signif- 
icance on the variable s in the range permitted by the on peak region. In this 
case we also plot the asymmetry and significance of the MSSM to compare 
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Silts 




Figure 7.6: On peak (M^/ — 3Tz' < M^+e- < Mz' + SF^/) asymmetry 
(left) and significance (right) versus the cut on rapidity Ycut- The anomalous 
charges are fixed to the general values Ql = 1, Qh^ = 1/2 and Qq = 3/4. 
The peak for the significance is for Ycut = 0.4. The position of the peak does 
not depend on the choice of the values for the charges, while the intensity of 
asymmetry and significance depends on it. 




Figure 7.7: Asymmetry (left) and significance (right) in MiAUSSM (blue) 
and MSSM (purple) around the peak of the Z' versus the invariant mass 
Me+e-- Vcut is set to 0.4, while the anomalous charges are fixed to the general 
values Ql = 1, Qh^ = 1/2 and Qq = 3/4 

the results. It is clear that the minimum in the asymmetry as function of s 
is due to the presence of the extra bosonic resonance, Z'. 

One side asymmetry and significance with fixed charges: the search 
of the peak for the significance 

The figure 17.81 shows the dependence of the on peak asymmetry and the sig- 
nificance versus the cut Ycut- It is clear that the significance has a maximum 
around Pz,cut = 580 GeV. Using this result we have obtained the figure 17. 9[ 
in which we keep fixed the cut and study the dependence of the asymmetry 
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Figure 7.8: On peak {Mz' - 3Tz' < M^+e- < Mz' + 3Tz') asymmetry (left) 
and significance (right) versus the cut on the longitudinal momentum Pz,cut- 
The anomalous charges are fixed to the general values Ql = 1, Qh^ = 1/2 
and Qq = 3/4. The peak for the significance is for Pz,cut = 580 GeV. The 
position of the peak does not depend on the choice of the values for the 
charges, while the intensity of asymmetry and significance depends on it. 




Figure 7.9: Asymmetry (left) and significance (right) in MiAUSSM around 
the peak of the Z'. Pz,cut is set to 580 GeV, while the anomalous charges are 
fixed to the general values Ql = I, Qhu = 1/2 and Qq = 3/4 

and significance on the variable s in the range permitted by the on peak 
region. 



Central asymmetry and significance with fixed charges: the search 
of the peak for the significance 

The figure [7!T0] shows the dependence of the on peak asymmetry and the sig- 
nificance versus the cut Y^ut- It is clear that the significance has a maximum 
around Ycut = 0.8. Using this result we have obtained the figure 17.111 in 
wich we keep fixed the cut and study the dependence of the asymmetry and 
significance on the variable s in the range permitted by the on peak region. 
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Figure 7.10: On peak (M^/ - 3Tz' < M^+^- < Mz' + 3Tz') asymmetry (left) 
and significance (right) versus the cut on the longitudinal momentum Pz,cut- 
The anomalous charges are fixed to the general values Ql = 1, Qhu = 1/2 
and Qq = 3/4. The peak for the significance is for Yc = 0.8. The position of 
the peak does not depend on the choice of the values for the charges, while 
the intensity of asymmetry and significance depends on it. 
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Figure 7.11: Asymmetry (left) and significance (right) in MiAUSSM around 
the peak of the Z'. Yc is set to 0.8, while the anomalous charges are fixed to 
the general values Ql = 1, Qh^ = 1/2 and Qq = 3/4 



Edge asymmetry with fixed charges: the search of the peak for the 
significance 

The figure 17.121 shows the dependence of the on peak asymmetry and the sig- 
nificance versus the cut Ycut- It is clear that the significance has a maximum 
around Ycut = 1-4. Using this result we have obtained the figure I7.13[ in 
wich we keep fixed the cut and study the dependence of the asymmetry and 
significance on the variable s in the range permitted by the on peak region. 
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Figure 7.12: On peak {Mz'-3Tz' < M^+e- < Mz' + 3Tz') asymmetry (left) 
and significance (right) versus the cut on the longitudinal momentum Pz,cut- 
The anomalous charges are fixed to the general values Ql = 1, Qh^ = 1/2 
and Qq = 3/4. The peak for the significance is for Yq = 1.4. The position of 
the peak does not depend on the choice of the values for the charges, while 
the intensity of asymmetry and significance depends on it. 




Figure 7.13: Asymmetry (left) and significance (right) in MiAUSSM around 
the peak of the Z'. Yq is set to 1.4, while the anomalous charges are fixed to 
the general values Ql = 1, Qh^ = 1/2 and Qq = 3/4 

7.4.2 Dependence of the asymmetry from couples of 
the charges 

We have already said that the asymmetry depends on the charges of the 
model. For the MiAUSSM we know that the free charges are Qq 
and Ql- In the previous section we have found the cuts that optimize each 
definition of asymmetry at the LHC Now, keeping them fixed to those values, 
we want to study the dependence of the asymmetry on the charges. 
Because we cannot represent on a single diagram a three variable function 
we chose to study each asymmetry with respect to the three possible couples 
of variables, i.e. Qh^, Qq, Qh^, Ql and Qq, Ql- 
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We have studied the value of the four asymmetries keeping alternatively one 
of the charges fixed to and varying the others two from —1 to 1. We limit 
the value in these ranges because in the SM all the charges are of this order. 
Furthermore we have found in sec 16.3. II that — 1 < Qh^ ^ 1; so for simplicity 
we have used the same region for the other two charges. 
In the following subsections we will show the contourplots of the results for 
each definition. 

RFB asymmetry with respect to the charges 




Figure 7.14: FB asymmetry (left) and significance (right) in MiAUSSM ver- 
sus Qh^ and Qq {Ql = and i/cut = 0.4) 
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Figure 7.15: FB asymmetry (left) and significance (right) in MiAUSSM ver- 
sus Qh^ and Ql{Qq = and i/cut = 0.4) 




-IJO -*J 'J.iJ Oj IJl -lit -OJ (1.(1 (tJ IJD 



Qq Qif 

Figure 7.16: FB asymmetry (left) and significance (right) in MiAUSSM ver- 
sus Qq and Ql{Qh^ = and i/cut = 0.4) 

In figures I7.14[ 17.151 and 17.161 we can see that Arfb is even for the ex- 
changes (Qh„, Qq, 0) {-Qh,,, -Qq, 0), {Qh^, 0, Ql) {-Qh^, 0, -Ql), 
(0, Qq, Ql) ^ (0, -Qq, Ql) and (0, Qq, Ql) ^ (0, Qq, -Ql). 
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One-side asymmetry with respect to the charges 




Figure 7.17: One-side asymmetry (left) and significance (right) in MiAUSSM 
versus Qh^ and Qq{Ql = and = 580 GeV) 
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Figure 7.18: One-side asymmetry (left) and significance (right) in MiAUSSM 
versus Qh^ and Ql{Qq = and p^^^^ = 580 GeV) 

In figures 17. 1 7[ l7.18l and l7.19W e can see that Aq is even for the exchanges 
(Qh„, Qq, 0) ^ {-Qh^, -Qq, 0), (Qh„, 0, Ql) ^ {-Qh^, 0, -Ql), 
(0, Qq, Ql) ^ (0, -Qq, Qi) and (0, Qq, Ql) ^ (0, Qq, -Q^). 
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Figure 7.19: One-side asymmetry (left) and significance (right) in MiAUSSM 
versus Qq and Ql{Qhu = and p^^^^ = 580 GeV) 

Central asymmetry with respect to the charges 




-IJO -Oj ilil (Jj IJl -1.0 -ilJ OjO Oj IJl 



Figure 7.20: Central asymmetry (left) and significance (right) in MiAUSSM 
versus Qh^ and Qq{Ql = and Yc = 0.8) 
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Figure 7.21: Central asymmetry (left) and significance (right) in MiAUSSM 
versus Qh^ and Ql{Qq = and Yc = 0.8) 




Figure 7.22: Central asymmetry (left) and significance (right) in MiAUSSM 
versus Qq and Ql{Qhu = Yc = 0.8) 

In figures 17.201 17.211 and 17.221 we can see that Ac is even for the exchanges 
Qq, 0) {-Qh^, -Qq, 0), (Qh^, 0, Ql) ^ {-Qh^, 0, -Ql), 
(0, Qq, Ql) -> (0, -Qq, Ql) and (0, Qq, Ql) -> (0, Qq, -Ql). 
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Edge asymmetry with respect to the charges 




Figure 7.23: Edge asymmetry (left) and significance (right) in MiAUSSM 
versus Qh^ and Qq{Ql = and Yc = 1.4) 
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Figure 7.24: Edge asymmetry (left) and significance (right) in MiAUSSM 
versus Qh^ and Ql{Qq = and Yc = 1.4) 
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Figure 7.25: Edge asymmetry (left) and significance (right) in MiAUSSM 
versus Qq and Ql{Qhu = and Yc = 1.4) 

In figures 17.231 17.241 and 17.251 we can see that Ae is even for the exchanges 
Qq, 0) ^ (-Qh„, -Qq, 0), (Qh„, 0, Ql) ^ (-Q/f„, 0, -Ql), 
(0, Qq, Ql) ^ (0, -Qq, Ql) and (0, Qq, Ql) ^ (0, Qq, -Ql). 
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7.4.3 Best fit for the Asymmetry in terms of the three 
free charges 

In the general case we do not impose any constraint except that the three 
charges can assume values between —1 and 1. From the cross section of our 
process, that can be found in section I7.3[ we can see that the amplitude is 
proportional to the fourth power of the charges. So, from the equations f l7.9p . 
(I7.25P and fl7.29p we know that the asymmetry must be a rational function 
in which both the numerator and denominator are fourth grade polynomials 
in the charges. 

^ El,,,=oa^AQHJiQQyiQL)' 
El,,,=ohAQHjiQQyiQL)' ^ ■ ' 

with i + j + k = n<4:. We have written a code that numerically calculates 
the coefficients of this fit. The simple idea is that by using the experimen- 
tal values for the on peak asymmetries and by considering only three of the 
four definitions we obtain a non-linear system with three equations and three 
variables Qq and Ql ) which could be solved numerically. In this way 

the asymmetry is useful for fixing the values of the U{iy charges. Moreover, 
once the values of the three charges are obtained by the previous system, 
the fourth definition of asymmetry can be used as a check for the validity of 
the model under exam. In fact its hypothetical experimental value must be 
recovered by using fl7.57p with the charge values already found, within the 
considered error (we use the medium relative error (MRE) for each asymme- 
try definition). As expected from the results of the previous section (although 
in that case one of the charges was fixed to 0), we have found out that the 
odd grade polynomials have negligible coefficients. In table [73] we have writ- 
ten the non zero coefficient of our polynomial for the four asymmetries. 

Note that this table contains only the statistical error and not the sys- 
tematic error due to the choice of the PDFs. 

We have calculated the goodness of the fit parameter (a perfect fit has 
= 1) and the medium relative error for these results, obtaining the re- 
sults that are showed in table 17. 4^ attesting the accuracy of the procedure. 
In particular the R^ value states (as we expect) that the errors in our fit are 
only due to numerical approximation in the calculation of the integrals that, 
given the cross section, give back the asymmetry . 
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Table 7.3: Coefficients of tfie fits for the four definitions of asymmetry 
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Ac 


Ae 




0,999 


0,999 


0,999 


0,999 


MRE 


0.008 


0.009 


0.019 


0.017 



Table 7.4: and Medium Relative Error for the polynomial fit of the asym- 
metry with respect to the three charges 
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Chapter 8 
Conclusion 



In this thesis we have studied many phenomenological aspect of the Mi- 
AUSSM. First of all, we have obtained the neutrahno mass matrix and we 
have derived the conditions to decouple the anomalous sector from the MSSM 
sector. Wc have found that at the tree level this condition is simply Qh^ = 0. 
In order to calculate the cross sections wc have calculated all the vertices of 
the model with and without mixing between the two sectors of the model. 
Then we have used these results to calculate the cross sections of the LSP- 
LSP annihilation and of the LSP-NLSP coannihilation in the case that the 
LSP is a Stiickelino and the NLSP is a MSSM neutralino. 
We have found that in the decoupled case the annihilation of two LSPs has 
a much weaker cross section than that needed to be in agreement with the 
WMAP results. So the Stiickehno cannot explain alone the DM relic density. 
Differently, we have found that if the Stiickehno coannihilates with an MSSM 
neutralino plus possibly other particles we can obtain the WMAP results. In 
the case of coannihilation with an NLSP that is almost an MSSM bino we 
have obtained that up to a mass gap of 10% we can satisfy the experimental 
constraints; in the case of coannihilation with an NLSP that is almost an 
MSSM wino and with its related chargino we find an agreement with the 
experimental results up to a 20% mass gap. 

We have also studied the general case, in which we do not decouple the two 
sectors of the neutralinos mass matrix. In this case we cannot have a pure 
Stiickelino as LSP, a pure bino as NLSP, etc. So the cross section of the LSP 
coannihilation is much more complicated that in the decoupled case. To deal 
with it we have modified the DarkSUSY package. 

Our motivation is simple: this package already contains routines to numer- 
ically calculate the relic density of particles in the most general cases. It is 
designed to work in an MSSM background, so we have added the new parti- 
cles, couplings and all the interaction vertices introduced by the MiAUSSM. 
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We have checked that our modifications do not affect the results given by 
the program if we shut down the extra U{1). 

We have used the modified version of DarkSUSY to extensively study the 
relic density of the LSP in the MiAUSSM with respect to the general set of 
free parameters given by the seven of the MSSM-7 and the five added by the 
anomalous extension. As in the case of decoupling of the anomalous and the 
MSSM sectors in the ncutralinos mass matrix, we have studied separately the 
case in which the NLSP is mostly a bino and the case in which it is mostly 
a wino. 

In both cases we have checked that the introduction of the extra U{1) does 
not lead to divergences or unphysical behaviour in the relic density We have 
also found that the relic density dependence on the MSSM parameters in 
our model is similar to that in the MSSM. For example we have found an 
example of funnel region, a phenomenon that is well known to happen in 
the MSSM if Maq ~ 2Mlsp and it is found in our model under the same 
condition. We have also checked that to satisfy experimental constraints we 
must obey to -1 < Qh^ < 1. 

Our main result in this context is the study of the dependence of the relic 
density on the anomalous masses Ms and Ma„- We have found that for mass 
gap of the order of 5% we satisfy the WMAP data in a region much wider 
than that in the case of mass gap around 10% and above. This implies that 
for small mass gap our model can satisfy the experimental constraints signif- 
icantly better than for bigger mass gap. 

In the second part of the thesis we have studied the asymmetry in the Mi- 
AUSSM for the process pp e~e~^ that is studied at the LHC. The asym- 
metry is used to impose constraints on theoretical models, to perform con- 
sistency checks and to calculate some quantities, usually the charges. 
We have calculated the elementary cross sections that contribute to this pro- 
cess, adding the contribution of the Z' to those of the 7 and the Zq, already 
present in the SM. Then we have calculated the asymmetry using four dif- 
ferent definitions that are viable at the LHC, namely rapidity dependent 
forward backward asymmetry, one side asymmetry, central asymmetry and 
edge asymmetry {A^ps, Aq, Aq and Ap respectively). 

All of these definitions contain a cut in the parameter space, necessary to 
obtain a result different from 0. We have calculated the value for the cut in 
each definition to have maximum significance and then we have studied the 
asymmetries (now "optimized " ) . 

The asymmetries and their related significancies are functions of the three 
free charges of the MiAUSSM. We have studied their dependence on the three 
possible couples of charges keeping each time one of them fixed. We have 
found that the asymmetries are even for exchange of sign of the charges as 
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we expected. 

We have also studied the general case, in which no charge is kept fixed. In 
this case we have calculated the value of the four definitions of asymmetry 
for discretized values of the charges in the interval [—1,1]. We have used 
this result to calculate a fit of the asymmetries as rational functions of the 
charges. Using these fits and using eventually founded experimental results, 
we have shown how can be obtained constrained values for the charges. 
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